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CHAPTER  1 


INTRODUCTION 

A  problem  of  great  interest  and  importance  in  underwater  acoustic 
signal  detection  is  the  coherent  combination  of  the  outputs  of  widely 
spaced  veceivers  to  form  a  very  large  array.  The  description  "widely 
spaced"  means  that  the  receivers  are  separated  by  distances  much  larger 
than  both  the  acoustic  wavelength  of  interest,  and  the  correlation  dis¬ 
tance  of  the  random  fluctuations  in  the  medium.  The  implication  is  that 
the  signals  received  by  the  individual  sensors  are  stochastically 
independent,  so  that  totally  new  methods  of  array  processing  are  re¬ 
quired  . 

There  are  two  aspects  of  this  problem,  which  involve  entirely 
separate  methods  cf  investigation.  The  first  is  pre-detection  coher¬ 
ent  combination,  by  which  the  receivers  use  a  priori  information  about 
the  state  of  the  medium  to  search  coherently  for  a  signal  source;  the 
objective  is  to  improve  signal  detectability,  and  performance  is  quanti¬ 
tatively  measured  by  the  array  gain.  Tne  second  aspect  is  post¬ 
detection  coherent  combination,  in  which  the  receivers  independently  de¬ 
tect  a  signal  source,  measure  the  signal  phase  in  real  time,  and  then 
form  a  coherently  focused  array  by  correctly  phase  shifting  the  signals. 
This  is  essentially  a  problem  in  signal  processing.  The  former  aspect, 
however,  is  primarily  a  problem  in  underwater  acoustic  wave  propagation 
in  a  random  medium,  and  requires  a  complete  analysis  of  the  space  and 
time  varying  character iotics  of  the  ocean  environment.  This  is  the 
problem  to  which  this  research  has  been  devoted,  and  which  is  the 


1 
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subject  of  this  dissertation. 


1.1  THE  PROBLEM  OF  COHERING  WIDELY  SPACED  RECEIVERS 


As  an  acoustic  propagation  medium,  the  ocean  presents  many  diffi- 


cult  problems  to  signal  reception.  The  speed  of  sound  underwater  varies 


in  space  and  time,  and  these  fluctuations  are  both  random  and  determin¬ 


istic.  Deterministic  spatial  variations  include  a  gradual  change  in 


the. sound  speed  with  depth,  causing  refraction  of  acoustic  rays  and  a 


multipath  signal  at  the  receiver.  In  addition-,  there  are  oceanographic 


phenomena  which  are  spac.e/time  random  processes,  resulting  in  unpredict¬ 


able  variations  in  the  sound  speed.  Some  of  these  fluctuations  are 


internal  waves,  tidal  phenomena,  currents,  eddies,  and  surface,  waves. 


The  combined  effect  of  these  sound  speed  changes  is  a  received  signal 


with  random  amplitude  and  phase  varying  spatially  and  temporally. 


Another  important  cause  of  spatial  phase  variations  is  the  effect  of 


ray  paths  which  change  with  range  and  the  resulting  spatial  change  of 


multipath  interference.  Additive  noise  further  degrades  signal  recep¬ 


tion;  the  primary  sources  are  ambient  noise,  which  is  random  and 


spatially  continuous,  and  discrete  noise  sources  such  as  shipping 


traffic  whose  characteristics  may  often  resemble  signals  of  interest. 


To  overcome  some  of  these  obstacles,  acoustic  sensors  are  combined 


into1  an  array.  An  amplification  and  a  phase  shift  are  applied  to  the 


received  output  of  each  sensor  and  the  results  are  summed.  If  each 


phase  shift  is  proportional  to  the  time  of  arrival  of  the  signal  at 


that  sensor,  then  the  array  is  phased  for  that'  particular  signal  source 


direction.  For  other  directions,  the  reception  will  be  partially 
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incoherent,  which  helps  in  rejection  of  noise.  In  a  conventional 
array,  the  spacing  between  sensors  is  on  the  order  of  a  wavelength. 
Since  the  correlation  distance  of  most  random  phase  and  amplitude 
fluctuations  is  much  greater  than  this,  each  sensor  sees  nearly 
identical  fluctuations  and  the  signal  outputs  of  the  sensors  can 
still  be  summed  coherently.  Therefore  the  primary  cause  of  degrada¬ 
tion  of  signal  reception  for  a  conventional  array  is  the  interfer¬ 
ing  noise.  Many  techniques  have  been  developed  for  noise  rejection 
and  can  be  found  in  the  literature. 

Another  important  function  of  an  array  is  localization  of  a  signal 
source.  A  measure  of  localization  ability  is  the  beamwidth,  which  is 
inversely  proportional  to  the  size  of  the  array  in  wavelengths.  The 
disadvantage  of  a  conventional  array  is  that  since  most  signal  sources 
are  in  the  far  field,  the  array  can  only  scan  in  angle;  range  informa¬ 
tion  must  be  estimated  from  the  intensity  of  the  received  signal. 

The  localization  is  then  limited  by  the  array  beamwidth.  However,  if 
the  receivers  are  separated  by  large' distances,  i.e.,  distances  very 
much  greater  than  a  wavelength,  and  on  the  same  order  of  magnitude  as 
the  range  of  interest  for  signal  detection,  then,  in  principle  some  of 
these  limitations  may  be  overcome.  The  array  could  then  scan  in  both 
range  and  angle,  since  signal  sources  would  be  in  the  near  field  of 
the  huge  aperture.  Also,  since  the  effective  beam  of  the  array  is 
then  a  very  small  two  dimensional  focal  spot,  resolution  ability 
would  be  greatly  enhanced. 

But  the  use  of  a  very  large  array  also  introduces  many  additional 
problems.  The  greatest  obstacle  is  that  of  localization  ambiguity. 


If  the  number  -of  component  sensors  is  small  and  they  have  omni¬ 
directional  reception,  then  there  are  numerous  locations  at  which  a 
signal  source  may  be  coherent  at  the  array,  and  localization  would 
be  impossible.  For  this  reason,  the  topic  considered  here  will  be 
limited  to  the  case  in  which  each  receiver  itself  is  an  array 
(henceforth,  reference  to  a  sensor  will  imply  a  subarray  receiver). 

This  limits  the  ambiguity  problem  to  the  area  of  overlap  of  the 
beams  of  the  subarrays,  before  coherent  combination.  Another  problem 
is  the  fact  that,  since  the  sensors  are  now  spaced  at  distances  much 
greater  than  the  correlation  lengths  of  randcm  oceanographic  fluctua¬ 
tions,  the  randomness  in  the  signal  is  independent  among  the  receivers. 
The  correct  phase  shift  to  apply  to  each  receiver  to  search  coherently 
for  a  signal  source  is  now  a  completely  unknown  quantity.  What,  if 
anything,  can  be  done  to  coherently  combine  these  receivers  to  form  a 
superarray  aperture,  and  thereby  improve  signal  detection  capability? 
This  is  the  question  which  will  be  addressed  and  answered  in  the 
dissertation. 

The  basic  approach  to  the  problem  is  as  follows.  A  beacon  signal 
source  is  placed  in  the  ocean,  and  radiates  a  known  waveform  to  each 
sensor.,  Each  sensor  measures  the  travel  time  of  the  signal  in  propaga¬ 
ting  through  the  random  ocean  channel.  This  information  yields  the 
correct  phase  shifts  for  the  superarray  to  focus  on  the  beacon.  The 
objective  is  to  scan  the  superarray  focal  spot  away  from  the  beacon  in 
search  of  a  signal.  But  since  the  ocean  is  fluctuating  both  spatially 
and  temporally,  the  distance  and  time  for  which  this  can  be  done  is 
limited  due  to  loss  of  coherence.  More  beacons,  will  be  required  so 
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that  the  superarray,  may  scan  from  beacon  to  beacon  to  maintain  an 
acceptable  level  of  coherence;  the  sensors  must  also  refocus  on  the 
same  beacon  as  often  as  is  determined  by  the  stability  time  of  the 
fluctuations.  The  results  to  be  presented  in  this  dissertation  will 
be  utilized  to  determine  these  required  beacon  spacings  and  refocus¬ 
ing  times,  for  specified  system  performance  parameters. 


1.2  THE  COHERENCE  FUNCTION  APPROACH  TO  THE  SOLUTION 

The  underwater  propagation  path  between  the  source  and  each 
sensor  is  modeled  as  a  random  channel  whose  stochastic  parameters 
depend  upon  the  oceanographic  fluctuations.  Since  the  random  variations 
in  the  received  signal  are  uncorrelated  among  all  sensors,  the  ability 
to  coherently  combine  the  distorted  signals  depends  upon  the  degree  of 
similarity  of  their  waveforms.  In  the  frequency  domain,  this  is 
viewed  as  a  measure  of  how  well  each  spectral  component  of  the  signal 
pairs  can  be  combined  in  phase,  despite  the  randomness. 

A  quantitative  measure  of  this  pairwise  coherence  is  given  by  the 
spectral  coherence  function.  Its  magnitude,  varying  between  zero  and 
unity,  is  the:  gain  in  received  signal  power  achieved  by’ combining  a 
pair  of  random  signals  with  partial  coherence;  a  value  of  unity  indi¬ 
cates  100%  gain  in  signal  power.  The  argument  of  the  coherence  function 
is  the  average  phase  difference  between  the  signals  necessary  to 
coherently  combine  them.  By  considering  all  possible  pairs  of  sensors 
in  an  array,,  the,  coherence  function  defines  an  important  array  perform¬ 
ance  parameter,  the  array  gain.  The  coherence  function  is  therefore 
the,  key  to  the  relationship  of  array  performance  to  oceanographic 
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fluctuations.  The  bulk  of  this  research  has  been  devoted  to  develop¬ 
ing  a  parametric  form  for  the  coherence  function,  which  can  be  used 
to  predict  array  system  performance.  The  general  theory  of  the 
coherence  function  is  presented  in  Section  2.4,  and  its  solution  for 
the  random  multipath  ccean  channel,  called  the  multipath  coherence 
function,  is  developed  in  Chapters  4  and  5. 

1.3  SUMMARY  OF  RESULTS 

By  means  of  the  model  of  uncorrelated  random  propagation  channels, 
an  expression  for  the  coherence  function  has  been  derived  in  terms  of 
the  parameters  of  real  oceanographic  fluctuations.  The  model  is 
generalized  to  include  scanning  distances  and  times.  Although  the 
results  include  the  most  recent  information  available  on  ocean  pheno¬ 
mena  such  as  internal  waves  and  tides,  the  structure  of  thev model  itself 
is  independent  of  these  data  and  can  easily  accommodate  future  changes 
or  new  theoretical  developments  in  oceanographic  fluctuations. 

The  results  of  the  analysis  demonstrate  a  simplification  that 
allows  numerical  results  to  be  computed  with  no  more  than  a  hand  calcula¬ 
tor.  The  derived  expression  for  the  multipath  coherence  function  is  a 
composite  of  three  factors  which  affect  signal  coherence:  determinis¬ 
tic  multipath  interference,  random  fluctuations  which  are  incoherent 
among  the  rays  of  a  multipath  set,  and  fluctuations  which  are  completely 
coherent  among  rays.  The  actual  multipath  configuration  can  be  ob¬ 
tained  from  a  ray  tracing  computer  program  or  from  experimental 
measurements  of  an  ocean  channel's  impulse  response.  The  second  fac¬ 
tor  is  dominated  by  internal  waves  and  the  spatial  variations  due  to 
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ray  paths  which  change  with  range.  The  last  factor  is  a  fluctua¬ 
tion  due  to  internal  tides.  The  value  of  this  mathematical  factori¬ 
zation  is.  that  it  permits  each  source  of  coherence  degradation  to  be 
analyzed  separately  and  the  relative  effects  of  each  to  be  compared. 

In  the  coherence-  function,  system  design  parameters  such  as  scan 

distance  and  scan  time  have  been  related  to  the  parameters  of  the 

ocean  fluctuations.  This  enables  a  determination  of  required  beacon 

spacings  and  beacon  refocusing  times  for  the  design  of  a  superarray 
* 

system.  These  results  are  then  applied  to  a  superarray  system  design 
to  demonstrate  practicality. 

Numerical  results  of  the  analysis  show  that  widely  spaced  receiv¬ 
ers  can  be  combined  with  partial  coherence  to  cover  large  ocean  areas, 
and  with  significant  realizable  array  gain.  In  addition,  it  is  shown 
that  such  a  system  design  is  practical  with  respect  to  required  density 
of  beacons  and  refocusing  times.  Methods  of  implementation  of  such  a 
system  are  proposed,  which  require  only  system  components  and  proce¬ 
dures  well  within  the  limits  of  current  capabilities,  both  technically 
and  economically. 

1.4  SUMMARY  OF  PREVIOUS  WORK 

The  primary  application  of  this  work  is  to  an  adaptive  array  tech¬ 
nique  known  as  self-cohering.  When  the  array  element  locations  are  not 
known-  accurately,  or  when  the  medium  has  a  randomly  varying  index  of 
refraction,  then  array  beamforming  and  scanning  must  be  performed  not 
by  a  priori  phasing  based  only  on  array  geometry  but  also  by  measurement 
of  signal  phase  from  a  direction  near  the  desired  source  location. 
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Self-cohering  techniques  applied  to  retrodirective  antenna  arrays 
were  first  discussed  in  [lj.  A  survey  of  current  and  previous  work 
in  self-cohering  techniques,  and  an  analysis  of  beatnforming  and  scan¬ 
ning  of  self-cohering  microwave  arrays,  is  given  by  Steinberg  [2], 

Most  of  the  current  research  in  self-cohering  techniques  for  very 
large  HF  and  microwave  arrays  is  being  done  at  the  Valley  Forge  Re¬ 
search  Center  [3], 

Self-cohering  techniques  for  arrays  were  first  introduced  into 
the  field  of  optics  in  the  early  1970's.  A  description  of  some  of 
this  work  can  be  found  in  [4].  In  principle,  the  techniques  are  identi¬ 
cal  to  those  used  for  antenna  arrays. 

Although  adaptive  techniques  have  been  used  in  underwater  acoustic 
array  processing  for  some  years  (5],  a  common  assumption  has:  been  per¬ 
fect  signal  coherence  across  the  array  aperture-.  A  discussion  of 
signal  processing  for  very  large  arrays  can  be  found  in  £6] ;  however, 
the  unlikely  assumption  of  perfect  signal  coherence  is  also  made  in 
that  report.  ■  ' 

The  most  important  aspect  of  this  work  is  the  analysis  of  signal 
coherence  in  random  ocean  channels.  There  are  two  different  definitions 
of  coherence  in  common  use.  In  the  field  of  electromagnetic  wave 
propagation,  particularly  in  optics,  the  measure  of  coherence  most 
commonly  used  is  simply  the  normalized  time-domain  cross-correlation 
function.  A  thorough  theoretical  analysis  of  the  significance  and 
application  of  this  coherence  function  is  given  by  Bcran  and  Parrent 
17],  and  they  also  give  a  survey  of  previous  researches.  It  is  surprise 
ing,  however,  that  all  of  these  researchers  were  unaware  of  the  other 
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definition  of  coherence  until  its  rediscovery  by  Mandel  and  Wolf  [8] 
in  1976;  First  introduced  in  time  series  analysis  by  Wiener  [9]  in 
1930,  it  is  defined  as  the  cross-power  spectral  density  of  two  time 
functions,  normalized  by  their  auto-power  spectral  densities.  Probably 
the  best  description  and  explanation  of  the  physical  significance  of 
this  spectral  coherence  function  is  given  by  Koopmans  [10] ,  who  also 
presents  a  complete  history  of  its  development.  Other  analyses  of 
this  coherence  function  and  its  use  can  be  found  in  Bendat  and  Piersol 
HU,  and  Jenkins  and  Watts  [12.], 

The  spectral  coherence  function  is  the  measure  which  is  used  in 
this  work.  Its  advantage  is  that  it  gives  an  unambiguous  quantitative 
measure  of  the  ability,  at  each  frequency,  to  coherently  combine  ran¬ 
domly  distorted  signals.  A  good  discussion  of  the  difference  between 
the  two  measures  of  coherence,  and  the  advantages  of  the  spectral  co¬ 
herence  function,  is  given  by  Roth  [13].  In  the  field  of  underwater 
acoustic  array  processing,  both  definitions  of  coherence  have  been 
used.  Use  of  the  cross-correlation  coefficient  in  the  definition  of 
array  gain  was  demonstrated  in  [14] .  Some  applications  of  the  spectral 
coherence  function  to  underwater  acoustic  processing  are  given  in  [15]. 

There  has  been  a  number  of  studies,  both  experimental  and  theoreti¬ 
cal,  of  coherence  of  acoustic  signals  in  a  random  ocean  environment 
and  its  effect  on  array  performance.  Smith  [16]  has  presented  an 
analysis  of  spatial  coherence  in  random  multipath  channels  duo  to  the 
effects  of  variations  of  multipath  interference  with  range.  However 
his  results  are  limited  to  separations  for  which  the  received  signal 
is  a  plane  wave,  and  the  random  variaii&ns  are  completely  correlated. 
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Jobst  and  Zabalgogeazcoa  [17,  18]  have  analyzed  the  effects  of  a 
moving  source  on  signal  coherence  in  a  multipath  channel.  Here  again 
the  signal  Is  assumed  to  be  a  plane  wave  across  the  array  and  the  phase 
fluctuations  are  also  assumed  to  be  completely  correlated  among  sensors. 
Hunk  et  al  [19]  have  determined  limits  on  coherent  processing  due  to 
phase  fluctuations  caused  by  internal  waves.  Their  analysis  is  also 
limited  to  small  sensor  separations  and  large  phase  fluctuations. 

The  major  difference  between  all  previous  work  and  the  work  to  be  • 

«  *  •> 

performed  here  is  that  the  former  has  been  limited  to  sensor  separations 
that  are  within  the  correlation  distance  of  the  random  fluctuations. 
Degradation  of  coherence,  then,  essentially  becomes  just  a  matter  of 
lack  of  correlation  between  the  randomness  in  signals.  But  this  gives 
no  insight  into  the  ability  to  combine  signals  with  partial  coherence 
when  the  receivers  are  far  beyond  this  correlation  distance.  If  the 
random  fluctuations  in  signals  received  by  widely  separated  sensors  are 
small  enough,  then  the  possibility  exists  for  achieving  some  gain  by 
properly  phase-shifting  one  signal  with  respect  to  the  other.  Another 
difference  from  previous  work  is  that  plane  wave  phase  shifts  are  gener¬ 
ally  used  for  conventional  bearaforming  and  scanning.  However,  these  do 
not  take  into  account  the  phase  bias  due  to  .multipath  and  oceanographic 
fluctuations.  By  using  the  true  average  phase  difference  between  sig¬ 
nals  as  predicted  by  the  coherence  function  in  terms  of  oceanographic 
fluctuations  a  further  increase  in  gain  may  be  realized.  The  spectral 
coherence  function  is  a  suitable  measure  of  this  potential,  and  it  is 
toward  this  end  that  most  of  this  research  has  been  directed. 


1.5  ORGANIZATION  AND  CONTENTS 

The  chapters  of  this  dissertation  are  organized  into  five  inter¬ 
related  levels  of  material  as  indicated  in  Fig.  1.1.  The  first  level 
consists  of  this  introductory  chapter  which  lays  the  groundwork  for  the 
dissertation  by  stating  the  problem,  the  approach  to  the  solution,  and 
giving  a  summary  of  results  and  previous  work.  The  second  level  is  com¬ 
posed  of  Chapters  2  and  3  and  presents  essential  background  information. 
Chapter  2,  "Underwater  Acoustic  Propagation  and  Array  Processing",  dis¬ 
cusses  the  wave  equation  and  ray  solution,  and  variations  in  the  sound 
speed  as  causes  of  phase  and  amplitude  fluctuations.  Some  characteris¬ 
tics  of  underwater  acoustic  signals  and  noise  are  presented  and  the 
general  theory  of  the  coherence  function  is  developed.  Basic  array  pro¬ 
cessing  theory  in  the  space  and  time  domains  is  discussed  including  the 
general  effect  of  the  randomness  of  the  medium.  Array  gain  and  its 
relationship  to  the  coherence  function  is  presented  and  methods  of  beam¬ 
forming  are  discussed.  Conventional  arrays  are  considered  with  respect 
to  their  characteristics  of  size,  directivity,  resolution,  and  correla¬ 
tion  of  random  fluctuations.  The  characteristics  of  very  large  arrays 
are  presented,  including  near  field  focusing  and  scanning,  resolution, 
uncorrelated  channels  and  uncorrelated  noise,  and  the  array  pattern. 

A  comparison  is  then  made  between  conventional  arrays  and  very  large  ar^ 
rays  (VLA) .  Finally,  the  topic  of  a  VLA  composed  of  conventional  sub¬ 
arrays  is  discussed. 

Chapter  3,  "Oceanographic  Fluctuations  and  Their  Effects  on  Propaga 
tion",  presents  the  characteristics  of  oceanographic  fluctuations  deter 
mined  from  experimental  observations,  and  classifies  them  according  to 
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CHAPTER  5  APPLICATIONS 


SUMMARY 


Fig.  1.1  Organization  of  dissertation. 
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their  effect  on  array  processing.  Types  of  fluctuations  are  then  dis¬ 
cussed  and  the  theory  of  those  fluctuations  which  are  relevant  to  the 
design  of  a  VLA  is  developed.  The  chapter  concludes  with  a  summary  of 
the  relative  importance  of  these  fluctuations  according  to  the  latest 
experimental  and  theoretical  results. 

The  third  level  of  organization  consists  of  the  central  theory  of 
the  dissertation,  presented  in  Chapter  4,  "The  Multipath  Coherence  Func¬ 
tion  for  Uncorrelated  Underwater  Channels".  This  level  makes  the  tran¬ 
sition  from  background  material  to  the  subject  of  the  thesis  and,  with 
few  exceptions,  follows  directly  from  the  first  chapter  for  one  thor¬ 
oughly  familiar  with  the  background  presented  in  Chapters  2  and  3. 

The  MCF  is  presented  as  a  new  measure  of  array  performance,  end  its 
physical  significance  is  explained.  The  MCF  is  derived  using  the 
stochastic  time-varying  channel  representation  of  multipath  propaga¬ 
tion  for-  general  oceanographic  fluctuations.  The  theory  is  then  extend¬ 
ed  to  include  the  effects  of  VLA  scanning  in  space  and  time.  The  re¬ 
sults  of  the  analysis  are  discussed  in  detail,  and  the  summary  presents 
a  prelude  to  the  development  of  the  MCF  in  terms  of  real  oceanographic 
fluctuations  in  Chapter  5. 

The  fourth  level  of  organization,  composed  of  Chapters  5  and  6,  is 
an  application  of  the  central  theory  of  Chapter  4  to  the  background 
material  presented  in  Chapters  3  and  2,  respectively.  Chapter  5, 

"The  Coherence  Function  in  Terms  of  the  Oceanographic  Fluctuations", 
incorporates  the  parameters  of  the  predominant  fluctuations  into  the 
MCF  and  analyzes  the  effects  of  each  on  coherence.  In  particular,  new 
theories  are  developed  for  the  effects  of  spatially  varying  multipath 
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interference  and  internal  tides.  Complete  numerical  results  are  given 
which  show  the  effects  of  source  range,  frequency,  multipath  character¬ 
istics,  and  scanning  on  coherence,  due  to  each  individual  source  of 
fluctuation.  Physical  interpretations  of  the  results  are  also  given. 

Chapter  6,  ’’Application  to  a  Super array  System  Design",  is  con¬ 
cerned  with  a  practical  application  of  the  previous  developments  to  the 
design  of  a  large  underwater  aperture  of  coherently  combined  subarrays, 
/h*  approach  to  a  complete  VLA  system  design  is  outlined,  including  such 
considerations  as  beacon  placement,  beacon  waveforms,  and  required  bea¬ 
con  spacings.  A  system  design  procedure  is  then  given  which  proposes 
a  methodology  for  implementation  of  system  specifications.  Finally, 
other  important  considerations  are  mentioned,  such  as  localization 
and  source  tracking. 

The  last  level  of  organization  is  comprised  of  Chapter  7,  "Summary 
and  Recommendations  for  Further  Study".  This  chapter  concludes  the 
work  with  an  interpretation  of  results  and  a  statement  of  all  limita¬ 
tions.  Recommendations  are  then  made  for  future  studies  of  relevant 
topics  not  considered  here. 

As  illustrated  in  Fig.  1.1,  the  essence  of  this  dissertation  can 
be  obtained  from  Chapters  1,  A,  and  7  which  contain  a  statement  of  the 
problem,  the  method  of  solution,  and  results,  respectively.  Chapters  2 
and  3  provide  a  basis  for  the  development  of  Chapters  5  and  6,  which 
latter  are  necessary  for  a  full  understanding  of  how  the  conclusions  of 
Chapter  7  follow  from  the  theory  developed  in  Chapter  A. 
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CHAPTER  2 


UNDERWATER  ACOUSTIC  PROPAGATION  AND  ARRAY  PROCESSING 

2,1  THE  WAVE  EQUATION  AND  RAY  SOLUTION 

The  propagation  of  an  underwater  acoustic  wave  obeys  the  wave 
equation  for  the  pressure 
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primarily  by  variations  in  temperature  and  pressure.  The  sound  speed 
Increases  as  temperature  and  pressure  increase,  resulting  in  a  sound 
speed  profile  as  shown  in  Fig.  2.1.  The  sound  speed  usually  has  its 
maximum  at  the  surface  where  the  temperature  is  the  highest.  The 
sound  speed  decreases  as  depth  increases  due  to  the  decreasing  tempera¬ 
ture  until  the  effect  of  increasing  pressure  causes  it  to  again  increase. 
The  depth  of  the  minimum  sound  speed  is  known  as  the  sound  channel 
axis.  Maximum  variations  of  c(z)  are  from  about  1480  m/sec  to 
1550  m/sec  and  depend  on  climate,  season  and  time  of  day. 

The  Eikonal  equation  is  valid  if 

X  Ag  «  1  (2.5) 

c 

i.e.,  if  the  fractional  change  in  the  sound  speed  gradient',  g  *  dc/dz, 
over  the  distance  of  a  wavelength  is  very  small  compared  to  f  *  c/X. 

The  surfaces  <{>(x,y,z)  =  constant  define  the  wavefronts  and  the  ray 
paths  perpendicular  to  these  wavefronts  can  be  found  once  c(z)  has  been 
specified.  An  example  of  ray  tracing  for  a  specific  sound  speed  profile 
is  given  in  Fig.  2.2.  From  the  Eikonal  equation  comes  the  underwater 
acoustic  equivalent  of  Snell's  law,  written  as 

c(z)  =  c(zv)cos9  (2.6) 

in  terms  of  the  sound  speeds  at  a  depth  z  and  at  the  vertex  depth  z , 
and  the  angle  0  which  a  ray  makes  with  the  horizontal  at  a  depth  z.‘ 

When  the  gradient  g  is  positive,  a  ray  is  concave  upward,  and  when  g  is 
negative,  it  is  concave  downward.  For  a  ray  which  leaves  a  source  at 


Fig.  2.2  Ray  paths  and  sound  channel 
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depth  z^  at  an  initial  angle  0q,  the  sound  speed  at  the  vertex  can  be 
found  from  (2.6).  The  ray  is  within  a  sound  channel  when  it  has  both 
an  upper  and  a  lower  vertex  and  all  rays  which  leave  the  source  at 
angles  smaller  than  0q  will  stay  within  this  sound  channel.  Certain 
rays  which  leave  the  source  will  reach  a  receiver  at  a  depth  z0  The 
underwater  sound  channel  is  therefore  characterized  by  multipath  propa¬ 
gation  between  source  and  receiver. 

The  total  phase  of  a  ray  in  propagating  from  source  to  receiver 
is  “  wT,  and  total  travel  time  is  found  directly  from  (2.4)  as 

T  “  \4fy  ■  (2-7> 

j 

ray  path 

The  travel  times,  pressure  amplitudes,  and  arrival  angles-  of  all  rays 
which  reach  a  receiver  are  usually  obtained  from  a  ray  tracing 
computer  program.  For  a  specified  sound  speed  profile,  source  range, 
frequency,  source  depth  and  receiver  depth,  the  program  will  compute 
the  above  quantities  for  all  possible  ray  paths  between  source  and  re¬ 
ceiver  . 

An  example  of  the  characteristics  of  multipath  propagation  is 
shown  in  Fig.  2.3.  The  source  and  receiver  are  separated  by  a  range 
R  =  500  km.  The  figure  shows  the  travel  times  and  relative  amplitudes 
of  the  rays  reaching  the  receiver,  all  of  which  are  bottom  reflected. 
The  nominal  or  average  travel  time  for  the  channel  is  seen  to  be  on 
the  order  of  T^  =  R/c  =  500/1.5  =  333  sec  and  is  called  the  bulk  time 
delay.  The  rays  arrive  in  pairs  with  approximately  the  same  amplitude, 


RELATIVE  AMPLITUDE 


Fig.  2.3  Ray  travel  times,  pressure  amplitudes,  and  arrival  angles 
(from  reference  9,  Chapter  3). 
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one  ray,  arriving  at  an  upward  angle  .and 'the;  other  at  the  same  angle 
in  a  downward,  direction.  The  initial  angles  of  these- rays  at  the 
source  are  also. equal  about  the  horizontal.  Rays  which’  leave  the 
'source  at  small  angles  haye  larger  amplitudes  due  to  shorter  path 
lengths  and  fewer  bottom  reflections.  As  the  initial  angle  of  each 
ray  pair  increases,  its  path  length  and  number  of  bottom  reflections 
increase,  and  the  relative  amplitude  decreases.  Due  to  the  attenua-r 
t'ion  of  high  angle  rays  over  long  ranges,  ray  arrivals  with,  signifi¬ 
cant  amplitudes  are  usually  limited  to  small  arrival  angles.  The  time 
between  the  first  ray  arrival  and  the  last  ray  arrival  is  called  the 
time  spread  of  the  channel,  Tg.  The  time  spread  usually  increases 
with  increasing  range  and,  for  the  figure  shown,  Tg  ■=  7  sec. 

The  received  pressure  field  for  a  multipath  channel  is  the 
superposition  of  K  individual  ray  arrivals  given  by 

K  ■  .  .  K  -jwT. 

H(o»)  -  l  P.  -  ejWt  £  A.e  '  k  ,  (2.8) 

k°l  K  k=l 


This  field  exhibits  interference  among  the  component  rays,  resulting  in 
frequency  selective  fading.  Depending  upon  the  arrival  times  and 
amplitudes  of  the  rays,  the  received  field  will  demonstrate  construc¬ 


tive  or  destructive  interference  at  different  acoustic  frequencies  as 
depicted’  in  Fig.  2.4.  The  received  field  will  be  at  a  maximum  at 
frequencies  for  which  the  rays  are  all  in  phase,  while  for  other  fre¬ 
quencies  it  may  fade  due  to  total  destructive  interference.  Frequency 
selective  fading  demonstrates  the  importance  of  a  frequency  domain 
analysis  of  multipath  channels. 
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2.2  EFFECT  OF  RANDOM  FLUCTUATIONS 

Besides  the  deterministic  variation  of  sound  speed  with  depth, 
there  are  oceanographic  fluctuations  which  are  random  processes  in 
space  and  time  and  cause  fluctuations  in  the  ocean  temperature,  result¬ 
ing  in  random  fluctuations  in  the  sound  speed.  Among  these  fluctuations 
are  internal  waves,  which  are  predominant,  internal  tides,  currents, 
and  eddies.  The  sound  speed:  is  now  given  by 

c*  *=  c (2)  +  6c(x,y,z,t)  (2.9) 

where  6c/c  is  typically  on  the  order  of  10  \  In  the  presence  of 
this  random  sound  speed  fluctuation,  the  rays  will  be  slightly  per¬ 
turbed  from  their  deterministic  paths  as  shewn  in  Fig.  2,. 5. 

The  effect  of  the  fluctuations  must  be  found  by  solving  the  wave 
equation  using  the  sound  speed  given  by  (2.9).  The  method  of  solution 
depends  upon  the  acoustic  wavelength,  range  from  source  to  receiver, 
and  the  correlation  lengths  and  times  of  the  random  fluctuations.  In 
this  T-rork,  the  solution  for  the  pressure  in  the  presence  of  ray  pertur¬ 
bations  will  be  restricted  to  the  geometrical  optics  region  where 
diffraction  effects  are  negligible  so  that  amplitude  fluctuations  are 
much  smaller  than  the  phase  fluctuations.  The  conditions  which  must  be 
satisfied  for  this  solution  to  be  valid  are:  . 

1.  The  wavelength  is  much  smaller  than  the  smallest  correla¬ 
tion  length  of  the  fluctuations, 

X  «  lQ  .  (2.10) 


MULTIPATH  PRESSURE  FIELD,  H(w) 
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Fig.  2.4  Frequency  selective  multipath  interference. 
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Fig.  2.5  Perturbation  of  ray  paths. 
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The  travel  time  is  much  smaller  than  the  smallest  correla¬ 
tion  time  of  the  fluctuations, 

T  «  tq.  (2.11) 

Diffraction  effects  (and  therefore  ray  amplitude  fluctua¬ 
tions)  are  negligible,  which  requires  that  the  size  of  a 
Fresnel  zone  be  smaller  than  the  smallest  correlation 
length  of  the  fluctuations. 


(2.12) 


For  homogeneous,  isotropic  fluctuations,  the  condition 
is 


(2.13) 


For  inhomogeneous,  anisotropic  fluctuations  such  as  internal 
waves,  the  diffraction  parameter  A  is  obtained  by  an  average 
over  a  ray  path, 


•  r 

(2.14) 


For  internal  waves  [1] 


7 


A  =  (50  Hz/f) (R/300  km).  -(2.15) 


The  total  mean  square  phase  fluctuation-  for  an  ind\ /idual 
ray,  4^,  satisfies  * 


A$>  <  l. 


(2.16) 


For  Internal  waves  [1]  v 

4>2  »’  (f/50  Hz)2(R/300  km).  (2.17,) 

«■  •  i  “ 

The  region  A  £  1,  $  <  1  corresponds  to  the  Rytov  solution 
of  the  wave  equation  .(the  method  of  smooth  perturbations)  in  which^ 
amplitude  variations  are  no  longer  negligible.  '  The  combination  of 
this  region  and  the  geometrical  optics  regime  comprises  the  unsaturated 
region,  in  which  a  propagating  wave  can  still  be  represented  by  an 
amplitude  and  a  phase.  This  is  .no-  longer  true,  however,  in  the 
saturated  regions  in  which  there  are  very  strong  perturbations  in  the 
ray  paths.  A  diagram  of  these  regions  for  internal  wave  fluctuations 
is  given  in  Fig.  2.6. 

With  the  restriction  to  the  geometrical  optics  region,  the  total 
til 

phase  of  the  ray  is 


*n'“ 


r  -dS  .  _  t  fdS  .  o>  \( 

|  c(z)-t-6c(x,y,z,t)  c  “  cn 

I  J  0  J 


ds  - 


j 

ktfl  ray  path 


(2,.  18) 


The  received  random  multipath  field  now  becomes 


,  .  jwt  p 
H(w)  *»  eJ  I  e  e 


(2.19) 


k=l 


where  t^  is  the  travel  time  variation  caused  by  perturbations  in  the 
ray  path.  The  -mean  square  phase  fluctuation  is 
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2.3  SIGNAL.  AND  NOISE  CllARACiERISTICS 

The  definition  of  signals  and  noise  is  somewhat  subjective  in 
that  it  depends  upon  what  type  of  acoustic  reception  is  of  primary 
Interest,  and  which  others  cause  interference  in  the  attempt  to  detect 
it.  A  signal  may  be  a  partially  coherent  narrow  band  acoustic  wave 
such  as  a  discrete  frequency  line  from  a  surface  ship,  while  the  noise 
may  be  incoherent  and  broadband,  such  as  ambient  noise  arising  from  a 
superposition  of  numerous  long  range  sources.  On  the  other  hand,  a  sig¬ 
nal  might  be  a  broadband  random  source,  while  interfering  noise  could 
be  narrow  band  and  highly  coherent  such  as  from  surface  ships.  In  this 
study,  a  signal  is  defined  as  any  acoustic  wave,  either  random  or  deter¬ 
ministic,  narrow  band  or  broadband,  which  originates  at  a  single  point 
source,  and  therefore  is  partially  coherent  at  separated  sensors.  Also 
the  noise  will  be  limited  to  random  broadband  ambient  noise  which  is 
incoherent  at  separated  receivers. 

2.4  THE  COHERENCE  FUNCTION 

Consider  an  acoustic  point  source  radiating  a  waveform  s(t)  which 
has  a  spectrum  S(to)  .  Assume  that  the  wave  propagates  without  attenua¬ 
tion  along  single  paths  to  two  separated  receivers.  In  each  channel, 
the  signal  incurs  a  time  delay  equal  to  its  travel  time,  a  random 
travel  time  fluctuation,  and  an  additive  noise.  The  travel  time 
fluctuation  is  slowly  varying  compared  to  duration  time  of  the  signal. 
The  received  outputs  are  then  spectrum  analyzed  and  summed  as  depicted 
in  Fig.  2.7. 

*ln  this  discussion  waveforms  are  truncated  aL  some  finite  time. 

Fourier  transforms  are  taken  over  this  finite  time  interval. 


FREQUENCY,  f (HZ) 
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s(t)  ^  S(w)  POINT  SOURCE 


Fig.  2.7  Signal  processing  for  random  channels 


30 


The  ensemble  average  power  output  of  this  resulting 
array  is  proportional  to 

pt  “  <isi  (a,)+Ni  (w)+s2  (w)+n2  (w)  i 

-  (|s1W|2/>+<|s2(«))|2)>  +2Re<^S1(u))S2*(ojy> 

+  ^|N1<W)|^+^|N2(U))|^  +2Re^N1(w)N2*(0))^  . 

Assuming  equal  xioise  power 

Nq(w)  e  | Nx <co)  ] *  ^|n2(w)|2^  , 

and  equal  signal  power, 

yw)  -  <(is1(w)i2^  -  (\s2^)\ -  y  , 


then 


PT  »  2S0(aj)+2S0(w)ReYs(w)+2N0(w)+2N0(a))ReYN(w) 


The  quantities  YgO1-')  and  Y^(w)  are  the  signal  and  noise 
functions,  respective!}',  defined  as 


V“> 


<^|s1M|2)!i<(|s2’ito)|2^>!s 


r  ,  0<|y„(w)|<1 


and 


/ni<w)N,  (w)y> 

V°  -  /:  N  ,  ,,  « — .  0<|v„w|<l 
OjM  2>  ’<(  » 2  («)  7J  • 


two  element 


(2.23) 


(2.24) 


(2.25) 


(2.26) 

coherence 


(2.27) 


(2.28) 


r*i£^W*£Z?TttXr’!B2T%XZ^^  1  *'-.’■ 
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Thn  significance  of  these  functions  is  more  apparent  in  their 
relationship  to  array  gain.  The  array  gain  is  defined  as  the  sig¬ 
nal  to  noise  power  ratio  of  the  array  divided  by  the  signal  to  noise 
power  ratio  of  the  individual  receiver. 


G  ,fA 

'  VB0 


(2.29) 


or,  equivalently,  as  the  signal  power  gain  of  the  array  divided  by 
the  noise  power  gain  of  the  array 


c  -  SA/S0 

‘  VS0 


Since  for  the  two  receiver  array, 


(2.30) 


SA  ‘  2S0+2S0ReYS 


(2.31) 


i  "$al 


and 


then 


na  "  2V2nor°yn  ‘ 


l+ReYc 


H-HeyN  * 


(2.32) 


(2.33) 


It  can  be  seen  from  (2.33)  that  Y„(u>)  is  a  quantitative  measure  of 

O 

the  gain  in  average  signal  power  achieved  by  combining  a  pair  of  sen¬ 
sors  with  partial  coherence;  a  value*  of  unity  indicates  100%  gain  in 


*zxarxz^tPm3m*r&ux>mxm 
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signal  power. 

In  order  to  make  the  significance  of  the  coherence  function  more 

clear,  it  will  be  assumed  that  the  random  travel  time  fluctuations, 

and  t^t  are  temporally  stationary  Gaussian  random  processes  which 

2 

have  zero  mean,  variance  O  ,  and  spatial  correlation  coefficient  p 
with  largest  correlation  distance  Lq,  so  that  for  receivers  with 
separations  greater  than  LQ,  p  =  0.  From  (2.27),  the  coherence 
function  is 


y3(U)  - 


-uVci-p)  -«VV 


*  (2.34) 


It  can  now  be  clearly  seen  that  the  coherence  depends  upon  both  the 
correlation  of  the  fluctuations  and  their  size.  The  most  important 
conclusion  to  be  made  is  that  if  p  =  0,  the  coherence  YoO*1)  is  not 

necessarily  zero,  and  in  fact  can  attain  values  very  close  to  unity  if 

2 

P  is  small  enough.  The  major  premise  of  this  dissertation  is  that  the 
random  fluctuations  are  stochastically  independent  due  to  the  large 
receiver  separations,  so  that  the  major  effort  is  directed  toward 
determining  the  size  of  the  random  fluctuation's.  The  above  expression 
also  hints  at  the  fact  that  the  argument  of  the  coherence  function  is 
the  average  phase  difference  between  the  signals  necessary  to  coher¬ 
ently  combine  them. 

For  arbitrary  signals,  the  coherence  function  is  formally  defined 


as  [3,  4,  5) 
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■Ys<“> 


G12<W) 


^(«)G2(«) 


(2.35) 


where  G-^Cw)  is  the  cross  power  spectral  density  of  the  received  sig¬ 
nals,  and  (w) ,  G^Cw)  are  the  auto  power  spectral  densities.  Its  two 
most  important  properties  are 

1  -  its  magnitude,  varying  between  zero  and  unity,  is  a  quantita¬ 
tive  measure  of  .the  ability  to  combine  random  signals  by  giving  the 
gain  in  average  signal  power. 

2  -  its  argument  is  the  average  phase  difference  between  the 
signals  necessary  to  coherently  combine  them  [4], 

It  should  be  noted  that  the  coherence  function  is  not  simply  the  fre¬ 
quency  domain  analog  of  the  normalized  tine  cross-correlation  function. 
The  correlation  function  is  normalized  only  to  the  mean  of  the  total 
power  in  each  channel,  but  the  coherence  function  is  normalized  at 
each  frequency  separately  [6].  Another  major  difference  is  that  the 
correlation  function  includes  the  entire  spectrum  of  frequencies  pre¬ 
sent  in  the  signal  waveform;  there  may  be  a  high  degree  of  coherence 
at  certain  discrete  frequencies,  but  this  information  will  be  lost  if 
coherence  is  low  over  the  major  portion  of  the  signal  spectrum.  This 
again  demonstrates  the  importance  of  frequency  domain  analysis. 


2.5  ARRAY  PROCESSING 

It  was  stated  in  the  previous  section  that  the  argument  of  the 
signal  coherence  function  is  the  average  phase  difference  between 


the  two  received  signals.  In  practice,  then,  each  sensor  pair 
would  phase  shift  the  received  signal  by  this  amount  before  adding 
the  receiver  outputs.  The  degradation  of  coherence  would  then  be 
determined  by  the  magnitude  of  the  coherence  function  which  is  a 
measure  of  the  random  phase  fluctuation  about  the  average. 

The  generalization  of  array  gain  to  an  array  of  N  sensors  with 
amplifications  and  phase  shifts  applied  to  their  received  signals 
before  combination  is  [7] 


N 

N  * 

l 

\  w  w 
m  n 

W“> 

m=l 

n=l 

N 

N  * 

l 

m=l 

)  w  w 
**,  m  n 
n=l 

W60 

(2.36) 

«v 

in  which  the  w^  are  the  complex  weights  for  the  amplifications  and 
phase  shifts,  and  the  signal  coherence  between  receivers  ra  and  n  is 

YSmn(w)»  With  YSnm(aj)  *  YSmn(w)  and  YSmm(w)  =  and  likewise  for 
YNmn^>*  Conventional  beamforming  is  the  choice  of  the  wffl  to  co¬ 
phase  for  the  average  signal  phase  difference  given  by  Ygmn(w). 

This  does  not  consider  the  effects  of  the  noise  on  array  gain. 

Adaptive  beaniforming  [3]  however,  consists  of  choosing  the  weight 
matrix  W  to  optimize  the  quantity  (2.36),  which  does  take  the  noise 

coherence  matrix,  I\.,  into-  account.  For  the  case  of  incoherent 

h 

noise,  becomes  the  unit  matrix,  and  the  two- methods  are  then  equiva¬ 
lent. 


Assume  that  the  sigr.al  coherence  has  the  same  magnitude  be¬ 
tween  all  pairs  of  sensors  so  that  |YSmnl  =  Yg  for  m  f  n,  and  that 
the  weights  have  unit  magnitudes,  with  their  phases  chosen  to  co¬ 
phase  perfectly  for  the  average  phase  difference  between  each  pair 
of  signals  given  by  the  argument  of  Ygmn*  If  the  pairwise  noise  co¬ 
herences  are  also  equal  so  lYjjnlnl  =  Y^  for  m  ^  n,  and  if  the  noise 
has  zero  mean,  then 


1+(K-1)YS 

1+(N-1)yn 


(2.37) 


This  expression  will  be  valuable  .in  comparing  conventional  arrays 
with  very  large  arrays. 

Another  important  performance  parameter  in  array  processing 
is  the  directional  power  response  of  the  array,  called  the  array 
pattern.  When  an  array  is  cophased  for  a  particular  source  direction, 
a  signal  arriving  from  a  different  direction  will  cause  a  different 
array  response  due  to  the  different  relati/e  path  lengths  among  sen¬ 
sors  to  the  new  source  direction.  The  direction  for  which  the  array 
is  cophased  is  the  primary  maximum  of  the  pattern  and  is  called  the 
main  beam.  For  some  other  directions  the  relative  path  lengths  will 
cause  a  partially  destructive  interference  resulting  in  a  region  of 
lower  power  response  called  the  sidelobe  region.  Besides  the  direction 
of  the  main  beam,  constructive  interference  will  occur  in  other 
directions  causing  additional  primary  maxima  in  the  pattern.  Since 
an  array  cannot  distinguish  which  primary  maximum  is  receiving  a 
signal  source,  this  results  in  source  location  ambiguities  sometimes 
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called  grating  lobes.  The  quantities  of  relevance  in  this  study 
are  the  main  beam  width  and  the  distance  between  primary  maxima. 

The  order  of  magnitude  of  these  quantities  can  be  estimated 
independent  of  noise  and  random  fluctuations. 

Consider  a  single  path  channel  with  constant,  non-random  sound 


speed,  and  a  source  at  a  range  R^.  The  signal  received  at  the  .nth 

sensor  is  proportional  to  s(t-T  )  which  has  a  spectrum  S  (w)  = 
-jioTn  n  n 

S(w)'-i  ,  where  T^  -  R^/c.  The  response  of  an  array  of  N  sensors 

to  this  signal  is 


N  H  *  * 

PA  c  l  l  w  w  S  (w)S  (oj) 
A  m=l  n=l  m  n  m  n 


*  “jw(T  -T  ) 

~  iS<«0|2IIwmw*e  »  n 


(2.38) 


To  cophase  for  a  signal  at  a  different  range,  R  the  complex 

jWT  nU 

weights  are  chosen  so  that  w  =  |w  |e  n  ,  with  T  „  =  R  „/c. 

n  n  nO  nO 

Eq.  (2.38)  then  gives  the  response  to  a  source  at  an  arbitrary 
range  R^;  when  R^  =  R^,  the  array  response  is  at  its  maximum. 


2.5.1.  CONVENTIONAL  ARRAY 

Consider  a  linear  array  of  length  L,  whose  N  receivers  are  located 
at  distances  d^  from  the  origin,  which  receives  a  signal  from  a  source 
at  range  R^,  as  depicted  in  Fig.  2.8.  A  conventional  array  is 
characterized  by 


1  -  receiver  separations  which  are  on  the  order  of  wavelengths, 
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(2.U) 


in  which  the  pattern  variable  in  sine  space  has  been  introduced. 


u  *=  sin0-sin8Q  . 


(2.45) 


The  width  of  the  main  beam  of  the  conventional  array  is  approximately 


Au_  -v  X 
B  L 


(2.46) 


and  the  spacing  between  primary  maxima  when  the  receivers  are  nearly 
equally  spaced  is  on  the  order  of 


4"m  '  &  • 


(2.47) 


Since  (2.45)  shows  that  the  maximum  range  of  u  is  2,  there  will  be 


no  ambiguities  if 


L  <  NX 
2 


(2.48) 


Eq.  (2.39)  implies  that  p  -  1  between  all  receiver  pairs,  so 
that  the  coherence  from  (2.34)  becomes 


rSmn  "  c 


-jk(d  ~d  )sin0 
J  v  m  n 


(2.49) 


which  shows  that  the  average  phase  difference  between  the  received 


signals  is 


<f>  .-(p  =  -k(d  -d  )sinO  . 

n  n  m  n 


After  cophaslng,  the  gain  of  (2.37)  for  unit  weights  is 


(2.50) 


y.”  "V  - 
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For  a  single  path,  non-random  medium,  =  R^/c,  and  the  array 
pattern  from  (2.38)  is 


.  jk(R  -R  ) 

p.  *=  n« w  e  n 

A  Llj  n  n 


(2.54) 


in  which  the  weights  should  be  chosen  to  cophase  for  the  desired 
source  location.  Unlike  the  conventional  array,  a  VLA  can  discrimin¬ 
ate  in  range.  It  is  therefore  convenient  to  give  the  pattern 

•k 

characteristics  in  units  of  length  in  both  range  and  azimuth.  The 
radial  width  of  the  main  beam,  called  the  depth  of  field,  is  on  the 
order  of 


(2.55) 


and  the  corresponding  width  of  the  main  beam  in  azimuth  or  cross¬ 
range  is 


AS_  ~  X  (Ro\  .  (2.56) 

\TI 

The  spacing  between  primary  maxima  when  the  receivers  are  nearly 
equally  spaced  is 


(2.57) 


in  range  and 

*Thcsc  results  were  obtained  from  computations  using  a  linear  VLA 
of  equally  spaced  receivers. 
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AS  .  ~  NXPO] 

M  vrJ 


(2.58) 


in  cross  range. 

It  is  apparent  that*  since  L  is  large  for  the  VLA  (on  the  order 
of  Rq) »  there  are  numerous  ambiguities  in  both  range  and  azimuth, 
with  spacings  on  the  order  of  wavelengths.  For  the  same  reason,  the 
VLA  beamwidth  is  much  smaller  tiian  that  of  a  conventional  array. 

The  characteristic  of  large  spacings  from  (2.52)  implies  that 
the  random  fluctuations  are  uncorrelated  between  receivers,  so  that 
pc  0.  From  the  simple  example  of  (2.34)  the  coherence  is  then 


22  -jw(T-T) 
-to  o  m  n 

e  e 


(2.59) 


The  coherence  now  depends  only  on  the  size  of  the  fluctuations 
determined  by  c f.  The  average  phase  difference  is 


-  -to(T-T)  . 
m  n  m  n 


(2.60) 


In  general,  this  phase  difference  cannot  be  predetermined  due  to 
the  random  fluctuations,  so  that  some  method  of  measurement  must  be 


used. 

Due  to  the  large  receiver  spacings,  the  ambient  noise  will  be 
incoherent,  so  that  =  0.  The  idealized  gain  from  (2.37)  then  be- 


C  »  l+(N-l)ys  . 


(2.61) 


comes 


For  Yg^O,  G->-I,  and  for  Yg‘,‘l»  G*Nt  its  maximum.  However,  in  con¬ 
trast  with  the  conventional  array,  there  is  now  no  limit  to  the 
attainable  gain  as  N  increases.  For  intermediate  values  of  Y<.,  and 
for  large  values  of  N,  G->NYg.  A  comparison  of  the  idealized  gain 
as  a  function  of  the  number  of  receivers,  N,  is  given  in  Fig.  2.10 
for  the  VLA  and  the  conventional  array. 

2.5. 2.1  VLA  OF  SUBARRAYS 

Consider  a  situation  in  which  there  are  N  indi/idual  omni¬ 
directional  receivers  with  which  to  design  an  array.  If  N  is  small 
then  it  is  not  practical  to  design  a  VLA  with  these  receivers  by 
separating  them  all  by  large  distances.  There  is  no  increase  in 
localization  due  to  directional  ambiguities,  and  gain  is  lost  due 
to  decrease  o'  signal  coherence  because  of  the  large  receiver 
spacings.  However  it  is  practical  to  subdivide  the  N  available 
receivers  into  coherently  combined  conventional  subarrays.  There 
will  be  an  increase  in  localization  ability  over  that  of  a  single 
conventional  array  of  N  receivers  since  each  subarray  has  a  beam 
which  can  intersect  those  of  the  other  subarrays,  and  the  ambigui¬ 
ties  of  the  VLA  are  limited  to  this  region  of  intersection. 

Consider  a  system  of  Ny  subarrays,  each  containing  Ng  receiv¬ 


ers.  The  subarray  gain  is 
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and  the  gain  of  a  VLA  of  omnidirectional  sensors  is 

Gv  -  1+<Kv-1)ys 

so  that  the  gain  of  a  VLA  of  subarrays  is 


For  yN  =  0, 


G  -  GyGs 


N, 


"  1+^NV“1)YS  *  l+(N_-l)Yk 


g  -  ns[i+(nv-i)ys1. 


(2.63) 


(2.64) 


(2.65) 


It  can  be  seer;  from  this  expression  that  even  small  values  of  VLA 
gain,  Gy,  can  be  very  significant.  For  example,  with  Ng  =20,  Ny  =  2, 
and  Yg  ~  .5,  from  (2.65),  G  =  30.  By  combining  only  2  subarrays  with 
a  coherence  of  only  50%  the  effective  number  of  elements  in  each  sub¬ 
array  when  they  are  used  incoherently  has  been  increased  from  20  to 
30.  The  expense  of  an  individual  subarray  system  including  its 
deployment,  operations  personnel,  signal  processing,  etc.,  may  be 
huge.  The  coherent  combinati.*  ,  oE  such  subarrays  requires  only  some 
additional  signal  processing  procedures  and  algorithms.  Therefore, 
from  a  cost  effectiveness  viewpoint,  a  VLA  gain  of  only  1.5  will  in¬ 
crease  the  value  of  such  a  large  system  by  this  same  factor,  with 


minimal  additional  expense. 

The  subarrays  may  still  be  used  incoherently  to  increase 


localization  ability  due  to  the  intersection  of  their  beams.  When 
they  are  combined  coherently,  the  localization  is  no  better  than  the 
incoherent  system,  but  the  value  of  the  increase  in  gain  achieved  may 
be  outstanding. 

Due  to  the  VLA  ambiguity  problem,  application  of  the  theory 
presented  here  will  be  limited  to  a  VLA  of  subarrays.  Further 
analysis  of  this  subject  will  be  presented  in  Chapter  6. 
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CHAPTER  3 


OCEANOGRAPHIC  FLUCTUATIONS  AND  THEIR  EFFECTS  ON  PROPAGATION 

3.1  CHARACTERISTICS  OF  OCEANOGRAPHIC  FLUCTUATIONS 

Oceanographic  fluctuations  cause  variations  in  the  amplitude 
and  phase  of  a  multipath  acoustic  signal.  Some  of  these  fluctuations 
are  environmental,  in  that  they  are  due  to  variations  in  tue  ocean 
medium  itself,  such  as  internal  wave  fluctuations,  independent  of 
the  presence  of  an  acoustic  signal. 

The  other  fluctuations  are  classified  as  acoustic,  since  they  de¬ 
pend  upon  the  presence  of  an  acoustic  signal  and  its  propagation 
characteristics.  Examples  of  this  type  are  spatial  multipath  varia¬ 
tions  due  to  changing  source  or  receiver  location,  and  frequency 
selective  fading  caused  by  multipath  arrivals  with  different  travel 
times.  In  addition,  the  environmental  fluctuations  cause  acoustic 
fluctuations,  since  signal  characteristics  are  influenced  by  the 
medium. 

Acoustic  fluctuations  may  be  spatial  and  temporal.  At  a  fixed 
location,  the  amplitude  and  phase  of  a  signal  will  vary  with  time, 
and  at  any  given  time,  they  will  vary  for  different  source  or  re¬ 
ceiver  locations.  A  good  example  of  environmentally  induced  acoustic 
phase  fluctuations,  which  demonstrates  their  spatial  and  temporal 
variability,  is  given  in  Fig.  3.1.  This  shows  the  results  of  a  13 

month  time  scries  of  acoustic  phase  taken  in  the  Straits  of  Florida, 
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and  reported  by  Steinberg,  et  al  [1J.  The  measurements  were  made 
at  three  fixed  colinear  hydrophones  at  ranges  of  7,  42,  and  43 
miles  from  a  fixed  source.  Besides  the  obvious  temporal  phase 
variation,  it  can  be  seen  how  the  phase  varies  with  receiver  separa¬ 
tion  at  a  given  time. 

Environmental  fluctuations  are  also  characterized  by  a  temporal 
spectrum  with  different  correlation  lengths  and  correlation  times. 
Periods  of  the  spectral  components  vary  from  minutes  to  months,  . 
and  characteristic  lengths  of  the  fluctuations  have  a  scale  ranging 
from  meters  to  thousands  of  kilometers.  It  can  be  stated  as  a 
general  rule,  that  the  lower  the  frequency  of  the  fluctuation,  the 
larger  are  its  energy  content,  correlation  length,  and  correlation 
time.  In  Fig.  3.1,  note  the  high  degree  of  correlation  between  the 
entire  time  series  of  H42  and  H43  due  to  their  small  separation, 
while  they  have  a  high  correlation  with  H7  only  at  the  longer 
period,  larger  amplitude  fluctuations. 

Environmental  fluctuations  can  also  be  classified  as  geographic 
and  non-geographic.  Non-geographic  fluctuations  are  those  which 
occur  in  all  areas  of  all  oceans  of  the  world,  such  as  interval  waves 
and  internal  tides.  Currents  and  eddies  are  examples  of  the  latter, 
and  occur  only  in  certain  areas  of  the  ocean  under  certain  conditions. 

An  excellent  report  on  current  knowledge  of  environmental  and 
acoustic  fluctuations  in  the  sea  and  measurement  techniques  is  pre¬ 
sented  by  Sykes  f 2 J .  This  report  summarizes  the  .results  of  measure¬ 
ments  done  over  the  last  decade  of  all  types  of  oceanographic 
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fluctuations  and  their  effects  on  propagation,  and  gives  a  com¬ 
plete  bibliography.  Any  reader  who  desires  further  information  con¬ 
cerning  oceanographic  fluctuations  should  consult  this  report. 

Further  analysis  here  will  be  limited  to  only  those  fluctuations 
which  are  relevant  to  this  study.  In  order  to  determine  this  limita¬ 
tion,  a  further  description  of  the  VLA  system  is  necessary. 

As  described  in  Section  1.1,  a  VLA  focuses  its  widely  spaced 
receivers  on  a  beacon  source  at  time  t.  Using' average  phase  shifts 
determined  by  the  pairwise  coherences  among  all  receivers,  the  VLA 
then  scans  a  distance  S  at  time  t  +  T.  The  statistics  of  coherence 
are  determined  by  considering  an  ensemble  of  identical  such  systems 
over  which  the  environmental  random  processes  of  interest  are 
stationary  in  space  and  time.  The  requirement  of  stationarity  first 
implies  that  each  member  of  the  ensemble  must  have  the  same  climate, 
meteorological  conditions,  and  season,  all  of  which  affect  the  nominal 
multipath  structure.  Secondly,  the  requirement  that 

T  «T0  (3.1) 

and 

S  «XQ  (3.2) 

where Tq  and  ^  are  characteristic  time  and  length  of  some  portion 
of  the  spectrum  in  time  and  space  of  all  environmental  fluctuations 
places  a  limit  on  the  fluctuations  which  must  be  considered  in  order  to 
maintain  stationarity.  For  T  on  the  order  of  hours  and  S  on  the  order 
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of  100  km  the  environmental  fluctuations  can  be  limited  to 

1  -  Internal  waves,  which  have  correlation  time  of  an  hour,  and 
correlation  length  of  several  kilometers. 

2  -  Internal  tides,  with  a  correlation  time  of  hours,  and 
correlation  distance  of  tens  of  kilometers. 

Larger  scale  fluctuations  with  on  the  order  of  days  or  longer,  and 
*o°f  hundreds  of  kilometers  or  larger,  can  then  be  omitted  and 
stationarity  will  still  be  maintained. 

In  order  to  maintain  a  uniformity  in  the  analysis  and  results, 
this  study  will  also  be  limited  to  those  types  of  fluctuations  which 
are  not  geographic  in  nature,  and  therefore  apply  to  all  oceans  of  the 
world.  The  analysis  thus  ignores  geographic  anomalies  such  as  currents 
and  eddies  which  may  further  degrade  coherence. 

3.2  TYPES  OF  OCEANOGRAPHIC  FLUCTUATIONS 

There  are  many  known  types  of  oceanographic  fluctuations,  and 
some  have  only  been  analyzed  experimentally.  Sykes  [2]  lists  the 
primary  causes  of  acoustical  fluctuations  as 

1.  Surface  waves  which  cause  frequency  spreading  of  the  signal 
spectrum  due  to  the  Doppler  effect.  Their  effect  is 
negligible  compared  to  other  fluctuations. 

2.  Internal  waves,  which  occur  due  to  varying  density  of  the 
ocean,  and  which  cause  variations  in  the  sound  speed.  They 
are  one  of  the  predominant  causes  of  acoustic  phase  fluctua¬ 
tions. 
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3.  Tidal  phenomena,  diurnal  and  semi-diurnal,  cause  changes 
in  water  depth  which  are  negligible  effects  for  deep 
ocean  propagation,  tidal  streaming  causing  currents  which 
are  a  geographic  effect,  and  internal  tides  which  are  one 
of  the  primary  causes  of  non-geographic  phase  fluctuations. 

4.  Rossby  waves  which  cause  long  term  large  space  scale 
fluctuations. 

5.  Solar  heating  which  causes  daily  changes  in  water  temperature 
and  acoustic  phase.  Its  effect  is  less  than  internal  tides. 

6.  Changes  in  lunar  declination  cause  large  phase  fluctuations 
with  a  period  of  27  days  and  a  large  space  scale. 

7.  Wind  influences  acoustic  phase  by  changing  the  water  tempera¬ 
ture. 

8.  Source  motion  causes  spatial  variations  in  multipath  inter¬ 
ference,  as  well  as  frequency  shifting  and  spreading  due  to 
a  different  Doppler  shift  for  each  ray  path. 

In  addition  to  these  from  reference  [2],  a  very  important  cause 
of  acoustic  phase  and  amplitude  fluctuations  is 

9.  Frequency  selective  fading  due  to  variations  in  multipath 
interference  as  frequency  varies.  This  effect  was  explained 
in  Section  2.1.  The  four  types  of  fluctuations  to  be  con¬ 
sidered  in  this  analysis  will  be  discussed  in  the  following 
sections  in  their  order  of  importance. 


I 
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3.2.1 


SPATIAL  VARIATIONS  DUE  TO  MULTIPATH  INTERFERENCE 


t 


As  source-to-receiver  range  varies,  the  travel  time  of  each 
ray  changes  at  a  different  rate.  This  causes  a  variation  of  the 
amplitude  and  phase  of  the  resultant  multipath  field  described  in 
Section  2.1.  For  large  changes  in  range,  the  number  and  types  of 
rays  which  reach  the  receiver  may  also  vary  due  to  changing  propaga¬ 
tion  geometry.  However,  for  smaller  range  variations,  the  ray  types 
and  number  of  arrivals  will  remain  constant.  This  latter  situation 


f 


i 


will  be  considered  here  for  simplicity;  in  any  event,  the  region  i 

| 

over  which  the  ray  characteristics  do  not  change  must  be  computed  jj 

from  a  ray  tracing  program.  | 

Clark,  et  al  [3],  have  analyzed,  through  a  ray  tracing  program,  ■ 

the  variations  in  resultant  phase  and  amplitude  for  a  source  moving  1 

from  500  km  to  520  km  at  various  speeds.  The  results  of  interest  to  | 


this  study  are  the  purely  spatial  variations  without  regard  to  the 
complicated  variations  due  to  the  Doppler  effect.  In  the  frequency 
domain  analysis,  the  effect  of  Doppler  shift  can  be  overcome  by  shift¬ 
ing  the  filter  frequency  of  the  receiver  by  the  proper  amount. 

Some  interesting  conclusions  can  be  drawn  from  the  results  of 
[3].  First,  there  is  a  linear  phase  trend  given  by  wT  where  T  =  R/c. 
When  this  effect  is  subtracted  out,  there  is  still  a  fluctuation  of 
the  resultant  amplitude  and  phase.  This  fluctuation  increases  as  the 
range  increases  from  the  reference  point.  Secondly,  as  the  reference 
range  increases,  the  spread  of  arrival  angles  generally  decreases, 
since  higher  angle  rays  are  attenuated  by  an  increasing  number  of 
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bottom  reflections.  This  implies  that  the  variation  in  resultant 
phase  will  be  less,  since  there  is  less  of  a  phase  difference  among 
rays  with  closely  spaced  arrival  angles. 

The  importance  of  these  variations  is  that  they  might  severely 
affect  scanning  ability  of  a  VLA,  since  average  phase  shifts  will 
be  used  to  scan,  and  there  may  be  large  variations  about  the  average 
due  to  the  spatial  multipath  interference.  Due  to  the  irapracticality 
of  computing  actual  variations  with  a  ray  tracing  program  for  each 
situation,  the  following  analysis  will  take  a  stochastic  approach 
to  the  solution. 


Theory 


Consider  the  expression  for  a  multipath  field  presented  in 


Section  2.1, 


K  ....  K  -jwT, 

H<w)  "  I  Pk  “  X  \e 
k»l  k=l  K 


(2.8) 


Each  of  the  K  rays  has  an  angle  of  arrival  6^.  Some  characteristics 
of  the  spread  of  angular  arrivals  are  symmetry  about  the  horizontal, 
and  a  rapidly  decreasing  density  of  arrivals  as  angle  increases  from 
the  horizontal.  If  (2.3)  represents  the  field  received 
from  a  source  at  range  R,  then  for  a  range  R  +  x,  where  x  is  small 
compared  to  R,  the  received  field  is  proportional  to 
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K  -juT. 

H(w)  =  Ue 

k=l 


X 

-jtojr  CO50j^ 


-  f  v'3“ 

k=l  * 


(3.3) 


The  quantity  t^  is  the  travel  time  variation  due  to  the  spatial 
changes  in  multipath  interference. 

It  is  desired  to  determine  how  the  amplitude  and  phase  of  H(w) 
vary  with  x  for  different  characteristics  of  the  arrival  angles,  0^. 
First,  assume  that  each  0^  is  an  independent  random  sample  from  some 
distribution  which  approximates  the  characteristics  of  the  determinis¬ 
tic  spread  of  0^.  Although  is  also  a  function  of  0^,  the  quantity 
of  interest  is  the  deviation  of  the  phase  and  amplitude  of  H(co)  from 
its  value  at  range  R,  regardless  of  the  values  of  T^,  so  that  the 
will  be  considered  to  be  non-random.  In  accordance  with  the  arrival 
angle  characteristics  stated  above,  the  ray  arrivals  will  be  approxi¬ 
mated  by  a  zero  mean  Gaussian  distribution.,  as  shown  in  Fig. 3. 2.  Since 
0^  is  small,  the  exponential  in  (3.3)  can  be  expanded  as 


-jojf  cos0k 
e 


0.  2 

.x  .  x  k 
“J“c  ~Y 


(3.4) 


The  expected  value  of  the  received  field  is  then 


K  -JojT. 

<H(w)>  =  Cc;(u))  l  A,  c  R  (3.5) 

k-l  K 
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where  Cg(<i))  is  the  characteristic  function  of  tg^> 

'  /  ”*^WtSkSs>  2  ^  1 

Cg(w)  *=  \e  /  =  [l+(2ao  )  ]  exp  j  -2ort-  -^tan 

(3.6) 

where  a  =  wx/2c  and  the  variance  is  O  =  n.  /  *  This  result 
shows  that  the  average  field  is  attenuated  as  its  resultant  phase  and 
amplitude  fluctuations  increase  due  to  increases  in  source  range  varia¬ 
tion,  frequency,  and  angular  ray  spread,  0.  In  addition,  the  result¬ 
ant  average  phase  is  a  composite  of  two  terms.  The  first  is  the 
nominal  phase  change  due  to  a  change  in  range  and  the  second  is  due 
to  the  spread  in  arrival  angles.  The  characteristic  function. 

Cg(w),  will  be  utilized  in  Chapter  5  to  determine  the  effects  of 
these  spatial  variations  on  coherence. 

3.2.2  INTERNAL  WAVES 

The  greatest  contribution  to  the  knowledge  of  internal  waves 
and  their  effect  on  acoustic  signals  has  been  made  by  oceanographers. 
Reference  [2]  gives  an  extensive  bibliography  concerning  work  in 
internal  waves. 

Internal  waves  arc  generated  in  regions  of  varying  density  in  the 


57 


ocean.  Propagation  of  the  waves  causes  random  variations  of  the 
density,  and  hence  the  sound  speed.  The  scale  sizes  of  internal 
wave  fluctuations  vary  from  meters  to  kilometers,  with  correlation 
distances  in  the  horizontal  much  greater  than  the  vertical,  i.e. 

Lg  »  Ly,  implying  that  the  ocean  is  anisotropic.  In  addition, 
the  sound  speed  fluctuations  caused  by  internal  waves  are  much 
greater  at  the  surface  than  at  greater  depths,  so  that  the  ocean  is 
also  inhomogeneous.  Internal  waves  are  also  characterized  by  a  dis¬ 
persive  spectrum;  roughly  speaking,  the  spectrum  of  the  phase 

—3 

fluctuations  varies  as  for  periods  ranging  from  1  hr.  to  24 

hr.  [4j . 

The  theory  of  internal  waves  used  here  will  be  based  largely 
on  references  [5]  and  [6],  This  theory  has  been  verified  by  compari¬ 
son  with  experiment  [4],  and  by  computer  simulation  [7].  Conclusions 
have  also  been  made  that  show  that  internal  waves  play  a  much  larger 
part  in  causing  acoustic  fluctuations  than  internal  tides  [8], 

There  are  three  important  quantities  which  characterize  the 
effects  of  internal  waves  on  acoustic  propagation: 

1,  The  strength  parameter,  <i>,  discussed  in  S'  ,  2.2,  which  is 
the  r.ra.s.  value  of  the  phase  fluctuation  for  a  single  ray 
in  the  geometrical  optics  region.  Depending  on  the  angle 
at  which  the  ray  crosses  the  sound  channel  axis,  it  has 
the  values  [5] 


(2.17) 
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axis  ray. 


(3.7) 


In  -order  to  make  the  frequency  dependence  explicit,  the  r.n.s. 
travel  time  fluctuation  is  introduced  as 


4>  “  4>/w 


(3.8) 


It  has  corresponding  values  given  by 


2  —8  2  —1 
4>  «*  (3.4x10  sec  km  )R,  steep  ray; 


(3.9) 


*=(6.8x10  ^sec^km  ^)R,  axis  ray.  (3.10) 


2,  The  diffraction  parameter,  A,  defined  in  Section  2.2. 

3.  The  phase  structure  function  defined  in  (2.21). 

For  a  horizontal  separation,  S,  at  constant  range,  R,  and 
a  temporal  separation,  T,  the  phase  structure  function 
for  internal  waves  is  [6] 


D(S,a)  «  2$‘ 


1 

/  s  \2  + 

1  f  T  \2 

2 

ui.4  kml 

2  1^1.6  hr j 

** 

(3.11) 


From  (2.21)  and  (3.1.1),  the  phase  correlation  coefficient 
for  internal  waves  can  be  deduced  as 


p(S,T)  =  1- 


S 


16.4  km  i 


+ 


1.6  hr 


(3.12) 
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The  structure  function  for  the  travel  time  fluctuationc 
is  then  written  as 

D(S,t)  =  D(S,T)/o)2 

-  24>2[l-p(S,T)3  .  (3.13) 

Internal  wave  fluctuations  are  such  that  they  cause  phase  fluctua¬ 
tions  which  are  uncorrelated  among  the  individual  rays  of  a  multipath 
field  [6],  Also  as  in  [6],  it  will  be  assumed  that  the  strength  para¬ 
meter  and  the  phase  structure  function  are  the  same  for  each  ray. 

3.2.3  -INTERNAL  TIDES 

Internal  tides  are  due  to  periodic  lunar  motion  and  cause  corres¬ 
ponding  periodic  variations  in  the  sound  speed.  There  are  two  pre¬ 
dominant  internal  tides,  the  semi-diurnal  and  the  diurnal.  In  the  deep 
ocean,  the  dominant  cause  of  tidally  induced  phase  fluctuations  is  the 
first  mode  M2  component  internal  tide,  which  has  a  period  of  12.42  hr 
and  a  wavelength  of  100  km.  The  internal  tide  propagates  outward  and 
inward  from  a  continental  shelf,  causing  a  sinusoidal  sound  speed  pertur¬ 
bation  with  the  same  wavelength  and  frequency  as  the  tide. 

An  acoustic  propagation  model  incorporated  in  a  ray  tracing  program 
by  Weinberg,  et  al  [9],  has  been  used  to  numerically  calculate  phase 
variations  d. e  to  internal  tides  based  upon  sound  speed  perturbations 
derived  by  Mooers  (10),  The  model  considers  an  acoustic  path  which  is 
perpendicular  to  the  direction  cf  propagation  of  the  internal  tide. 
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The  results  confirm  that  there  are  no  marked  differences  in  the  phase 
behavior  for  different  ray  paths,  and  that  phase  fluctuations  due  to 
internal  tides  can  therefore  be  considered  as  coherent  among  the 
individual  rays.  Since  the  phase  behavior  is  independent  of  the  ray 
considered,  it  is  sufficient  to  restrict  the  analysis  to  a  ray  on  the 
sound  channel  axis,  and  to  assume  that  it  yields  a  good  description  of 
the  bulk  time  delay  variations. 

Fig.  3.3  depicts  the  geometry  of  an  axis  ray  propagating  from  a 
range  R  at  an  angle  ({>  with  respect  to  the  wave  normal  of  the  internal 
tide.  The  axis  sound  speed  at  some  range  R  from  the  receiver  varies 
according  to  the  tidal  propagation  as 


c(r,t)  =  c0+Ac0sin(wTt-kTrcos(f>) 


(3.14) 


where  c^  is  the  unperturbed  axis  sound  speed,  Ac^  is  a  small  perturba¬ 
tion  due  to  the  internal  tide,  and 


WT  *=  2tt/ (12,42  hr.), 
kj  *=  2ir/ (100  km. ) 


(3.15) 

(3.16) 


are  the  radian  frequency  and  wavenumber,  respectively,  of  the  M2  tide. 
The  travel  time  of  the  ray  is  given  by 


T  - 


■R 


dr 


c(r,t) 


(3.17) 
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.2  Distribution  of  ray  arrival  angles. 


INTERNAL  TIDAL 


3  Acoustic  propagation  geometry  for  internal  tide 
fluctuations. 


Since  ACq  «  Cq,  the  result  of  integration  simplifies  to 


T  *=  T, 


,  slnfo..t-|t.RcoS»/2)sin(kXR,:osli’/21 
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(kTRcos«/2) 


where  T^  =  R/c^  *s  fc^e  travel  time  in  the  absence  of  the  internal  tide. 

Some  important  observations  concerning  the  travel  time  fluctua¬ 
tions  can  now  be  made,  based  upon  the  above  expression.  The  maximum 
variation  occurs  when  the  acoustic  path  is  perpendicular  to  the  direc¬ 
tion  of  the  internal  tide  propagation,  i.e.,  when  the  acoustic  signal 
propagates  parallel  to  a  continental  shelf;  the  minimum  variation  is 

when  the  acoustic  path  is  in  the  same  direction  as  the  internal  tide 

*  A 

(4>  =  0).  This  is  the  opposite  of  the  claim  made  in  reference  [9]  . 

Secondly,  it  can  be  seen  that,  for  very  long  source  ranges,  the 

fractional  variation  in  sound  speed  decreases. 

The  model  shows  excellent  agreement  with  experiment  [9]. 

3.2.4  FREQUENCY  SELECTIVE  MULTIPATH  INTERFERENCE 

The  interference  of  multipath  arrivals  with  different  travel  times 
causes  an  acoustic  fluctuation  in  the  frequency  domain  called  frequency 
selective  fading  which  was  briefly  described  in  Section  2,1.  This  is 
listed  as  the  least  important  acoustic  fluctuation  to  be  considered 


*There  were  no  computations  done  in  this  reference  for  the  case  of  acous¬ 
tic  propagation  in  the  same  direction  as  internal  tide  pr jpngation ,  which 
would  require  a  range  dependent  sound  speed  profile. 
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because  it  is  a  semi-periodic  function  of  frequency,  while  the  other 
fluctuations  increase  monotonically  with  frequency.  However,  it  plays 
an  important  role  in  the  analysis  of  coherent  frequencies  and  coherent 
bandwidths  which  will  be  discussed  in  Chapter  5. 


3.3  SUMMARY 

A  survey  and  classification  of  oceanographic  fluctuations  has  been 
presented.  In  order  to  maintain  a  uniformity  in  applications  of  the 
results,  consideration  of  environmental  fluctuations  has  been  limited 
to  those  v/hich  are  not  geographic.  However  this  does  not  preclude  the 
later  inclusion  of  anomalous  fluctuations,  since  the  multipath  coherence 
function  developed  in  Chapter  4  will  have  general  applicability  because 
of  a  classification  of  fluctuations  according  to  those  which  are  com¬ 
pletely  correlated  among  rays  (e.g.  internal  tides),  and  those  which  are 
uncorrelated  (e.g.  internal  waves  and  spatial  variations). 

The  justification  for  considering  only  internal  waves  and  tides 
as  the  predominant  types  of  environmental  fluctuations  is  due  to  the 
very  large  scale  sizes  and  correlation  times  of  other  fluctuations 
relative  to  VLA  scan  distances  and  scan  times.  In  principle  the  theory 
could  be  extended  to  larger  systems  which  must  consider  these  fluctua¬ 
tions  if  more  was  known  about  their  characteristics.  However  the  much 
larger  amplitude  of  these  fluctuations  would  make  the  design  of  a  larger 
system  impractical,  so  that  the  size  of  a  VLA  system  would  still  be 
determined  by  the  smaller  fluctuations  considered  here.  In  addition, 
the  combined  effect  of  the  smaller  fluctuations  on  coherence  is  large 
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enough  to  preclude  consideration  of  larger  fluctuations. 

There  has  been  some  controversy  between  the  oceanographic  and 
acoustic  communities  concerning  the  relative  importance  of  internal 
vaves  and  internal  tides.  A  paper  describing  a  recent  experiment 
claims  that  70%  of  the  energy  in  phase  fluctuations  of  periods  less 
than  one  day  is  due  to  the  semi-diurnal  internal  tide  [11].  However, 
the  large  frequency  bandwidth  used  in  making  that  conclusion  includes 
a  large  portion. of  energy  due  to  high  amplitude  internal  wave  fluctua¬ 
tions,  while  the  internal  tide  itself  has  an  extremely  narrow  bandwidth. 
An  analysis  using  uniformly  accepted  values  for  sound  speed  fluctua¬ 
tions  due  to  both  internal  waves  and  internal  tides  has  shown  that  90% 
of  the  total  energy  ir.  the  phase  fluctuations  is  due  to  internal  waves 
[8].  Internal  waves  therefore  have  the  larger  effect  on  phase  fluctua¬ 
tions  and  it  will  be  shown  in  Chapter  5  that  internal  tides  have  a 
negligible  effect  on  coherence  compared  to  internal  waves  and  spatial 
fluctuations. 
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.CHAPTER  4 


THE  MULTIPATH  COHERENCE  FUNCTION  FOR  UNCORRELATED  UNDERWATER  CHANNELS 

4.1  INTRODUCTION 

This  chapter  introduces  a  new  measure  for  determining  the 
coherence  of  acoustic  signals  in  multipath  channels  which  have  random 
fluctuations  that  are  uncorrelated  between  channels.  This  multipath 
coherence  function  (MCF)  is  based  upon  a  formulation  of  the  spectral 
coherence  function  in  terms  of  the  random  multipath  transfer  functions. 
The  MCF  allows  each  channel  to  be  analyzed  individually,  and  separates 
the  effects  of  random  fluctuations  from  the  effects  of  deterministic 
multipath  interference  (frequency  selective  fading). 

The  physical  significance  of  coherence  was  explained  in  Section 
2.4.  Coherence  is  a  quantitative  measure  of  the  extent  to  which  it  is 
possible  to  combine  randomly  distorted  signals  in  phase,  at  each  fre¬ 
quency  in  the  signal  spectrum.  The  coherence  is  quantitatively  related 
to  the  array  gain  in  that  it  is  a  measure  of  the  increase  in  received 
signal  power  achieved  by  combining  signals  with  partial  coherence 
relative  to  combining  them  incoherently  (i.e.,  adding  intensities). 

All  previous  analyses  of  coherence  have  been  limited  to  the 
situation  in  which  the  receivers  are  located  within  the  correlation  dis¬ 
tance  or  "patch  size"  of  the  random  fluctuations.  Most  of  these  inves¬ 
tigations  have  used  this  correlation  length  as  the  limiting  sensor 
separation  for  which  coherent  processing  can  be  performed.  Smith  [1] 

has  presented  an  analysis  of  spatial  coherence  in  random  mulLipath 
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channels  due  to  the  effects  of  spatial  variations  in  multipath 
interference.  However,  his  results  are  limited  to  separations  for 
which  the  received  signal  is  a  plane  wave,  and  he  assumes  that. random 
variations  are  large,  and  completely  correlated  between  sensors. 

Jobst  [2]  has  analyzed  the  effects  of  a  moving  source  on  signal  co¬ 
herence  in  a  multipath  channel  by  assuming  the  number  of  ray  arrivals 
to  be  a  random  variable.  Here  again,  the  signal  is  assumed  to  be  a 
plane  wave  across  the  array,  and  phase  fluctuations  are  assumed  to  be 
completely  correlated  between  sensors.  Munk,  et  al  [3]  have  deter¬ 
mined  limits  on  coherent  processing  due  to  phase  fluctuations  caused 
by  internal  waves  whose  characteristics  they. have  thoroughly  analyzed 
[4,  5].  Their  analysis  also  is  limited  to  snail  sensor  separations, 
and  their  criterion  for  degradation  of  coherence  is  not  quantitatively 
related  to  array  gain.  Beran  and  McCoy  [6,7]  have  done  analyses  of 
coherence  in  ocean  channels  using  the  mutual  coherence  function. 

Again  their  work  is  limited  to  plane  wave  propagation  within  the 
correlation  distance  of  the  fluctuations. 

There  are  two  major  differences  between  all  known  previous  work 
and  the  results  to  be  presented  here;  the  former  have  all  been  limit¬ 
ed  to  the  case  in  which  the  sensor  separations  are  small  enough  that 
they  are  within  the  correlation  distance  of  the  random  fluctuations, 
and  each  ray  defines  a  plane  wave  arrival  across  the  sensors.  The  re¬ 
sults  in  this  dissertation  apply  when  the  receivers  have  uncorrelnted 
fluctuations,  and  each  may  even  receive  an  entirely  different  multipath 
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4.2  DERIVATION  OF  THE  MULTIPATH  COHERENCE  FUNCTION 

Quantitatively,  the  coherence  function  can  be  defined  in  terms 

of  the  power  spectral  densities  of  the  received  signals,  using  the 

stochastic  time  varying  channel  approach  [8,  9,  10].  Consider  a 
"M 

point  source  radiating  a  signal  s(t)  with  spectrum  S (to)  which  propa¬ 
gates  through  two  linear,  random  multipath  channels  as  shown  in 

Fig.  4.1.  Since  the  channels  are  time  dispersive,  the  impulse  re- 

i* 

sponse  is  of  the  form 


K 


h(t)  =  ][  Ak6(t~Tk) 


(4.1) 


in  which  K  is  the  number  of  ray  arrivals,  Ak  is  the  amplitude  of  a 
ray,  and  Tk  is  its  travel  time  including  random  fluctuations  which  are 
slowly  varying  compared  to  signal  duration  time  and  travel  time. 

The  transfer  function  is  proportional  to 


K  -jwT 
H(to)  *  l  Ke 
k=l  K 

The  resulting  output  spectra  at  sensors  m  and  n  are 


(4.2) 


S  (to)  =  H  (u)S(to) 
m  m 


(4.3) 


and 

sn(w)  -  Hn(to)S(w), 


(4.4) 


*This  insures  that  the  source  is  coherent. 

fReflcction  phase  shifts  have  been  omitted.  They  will  only  affect  the 
exact  locations  of  coherent  frequencies  (Section  5.2.4)  which  must  be 
found  by  measurement. 
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in  which  H  (10)  and  H  (to)  may  be  different, 
m  n 

It  is  desired  to  coherently  combine  the  received  signals 
S^co)  and  Sn(io).  A  measure  of  the  ability  to  do  so  is  given  by  the 
spectral  coherence  defined  in  Section  2. A  as 


W“>  ■ 


G  (to) 
mn _ 

/G  (w)G  (to) 
m  n 


(2.35) 


Since  the  complex  transfer  functions  of  the  channels  are  random,  it  can 
easily  be  shown  that 


II 

3 

j 

/h  (io)H*(o>)' 
\  m  n  j 

)  G(to) 

(4.5) 

and 

G  (co)  = 
m 

<IV»I*> 

G(to) 

(4.6) 

where  { •  /  denotes  an  average  over  an  ensemble  of  random  processes 
as  described  in  Section  3.1,  and  G(co)  is  the  pov.’er  spectral  density 
of  the  input  signal,  s(t).  The  coherence  can  then  be  written  as 


Ww)  - 


<Hm(w)/(o>))  G(w) 
<(|Hm(w)|2^>H  <|h>)|2>*  G(w) 


(4.7) 


The  doherence  therefore  is  independent  of  the  input  signal  and  depends 
only  on  the  properties  of  the  channel.  If  the  random  transfer  functions 
of  the  channels  are  independent,  the  multipath  coherence  function  can 


be  written  as 
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Yn(B)Yn(M) 


(4.8) 


where  Y  (to)  and  y  (u)  will  be  called  the  auto-coherences, 
'm  n 


The  significance  of  this  result  is  very  important.  First,  it 
demonstrates  the  existence  of  partial  coherence  when  the  channels  are 
uncorrelated.  Second,  the  convenient  factorization  into  two  auto¬ 
coherences  allows  each  channel  to  be  analyzed  independently  of  all  the 
others.  This  implies  that,  for  an  array  of  N  receivers,  only  N  auto¬ 


coherences  must  be  computed  to  completely  determine  array  gain.  This 


,  N(N~1) 

can  represent  a  great  savings  compared  to  the  computation  of  — ^ - 


much  more  complicated  pairwise  coherences  if  the  channels  are  net 
independent.  Although  the  most  important  oceanographic  fluctuations, 
i.e.  internal  waves,  are  independent  among  receivers  of  a  VLA,  the 
MCF  can  easily  be  generalized  to  include  an  additional  type  of 
fluctuation  which  may  have  some  degree  of  correlation  between  channels 
e.g.  internal  tides.  The  effect  of  this  generalization  will  be  the 
addition  of  a' third  factor  to  the  MCF  which  is  the  coherence  due  to 
the  correlated  fluctuations  alone. 

The  random  travel  time  of  a  ray  will  now  be  written  in  terms- of. 
its  components  as 


Tk  "  wv  rny  k 


(4.9) 


identified  as: 
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Tk0  -  the  nominal  travel  time  of  the.  ray  in  the  absence 
of  any  fluctuations. 

t^  -  a  zero-mean  fluctuation  which  is  independent  and 
identically  distributed  among  the  rays  of  a  chan¬ 
nel  and  uncorrelated  between  channels, 
t^  -  a  fluctuation  which  is  completely  correlated  among 
rays  of  s  channel,  having  the  same  value  for  each 
ray.;  there  may  be  some  degree  of  correlation  between 
channels,  and  it  is  independent  of  the  fluctuation 


The  transfer  functions  of  the  two  channels  are  therefore 


H  <u>) 
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The  numerator  of  the  MCF  is  then 


<1lm(»)H*(U)>  - 


4  *  ‘tTn  tln>>  *  k,  (“>  <pw>u  (4  L  («» <Kn)>>n0  <u>  • 

k  k 

.-jw(t  -u  )V 

<e  1U  U1  X  (/ 


“jwTkOm  3wT£0n 


(4.12) 


where  (to)  and  H^qCu)  are  the  normalized  transfer  functions  in 

the  absence  of  fluctuations,  and  c,.  (w)  is  the  characteristic  func- 

wm 

.  tion  of  ttJ.  .  Similarly, 


kx, 


^aL+ 


CH»(U> 


(A. 13) 


with  an  analogous  expression  for  channel  n.  Tn i  ratio  of  coherent 

field  intensity  to  incoherent  field  intensity  is  the  quantity 
2  2 

(Ja^)  When  the  ray  amplitudes  are  equal,  this  ratio  is 

k  k  1011 

equal  to  K^,  the  number  of  rays  in  channel  m.  Henceforth,  the  para¬ 
meter  will  be  substituted,  with  the  understanding  that  it  designates 
this  ratio  when  the  amplitudes  are  unequal.  The  square  magnitude  of 
the  MCF  can  then  be  written  as 


.  (u)|2_  V^<“>lHmO<M>l2  . 

s“”  i+t*i.i”«o<“>i2-i)4(“)  i+t\iHn„(“)i2-u4n(“) 


kX  kl  A  \ 

m1  mO 


lYm(w)|  |Yn(w)| 


(A.1A) 


where  Y  (w),  y  (w)  are  the  auto-coherences.  It  is  shown  in  the 
ra  n 

appendix  that  each  of  these  factors  has  an  envelope  given  by 
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^-jw(tTm-tTn)^ 

yif 

•  W“)W">W">YMn(“)W“>  ' 


The  MCF  has  conveniently  factored  into  five  terms  that  permit  the 
effects  of  the  random  fluctuations  to  be  analyzed  independently  of 
the  effects  of  multipath  interference  as  can  be  seen  by  writing 


YSmn  =  YWYTYM 


(4.17) 


in  which  the  effect  of  uncorrelated  ray  fluctuations  is 

Y  e  Y  Y  > 

'W  Wm  Wo 

that  of  correlated  ray  fluctuations  is 
YT  =  YTmn  ’ 

and  the  effect  of  deterministic  multipath  interference  is 
YM  =  YMniYMn  * 


(4.18) 


(4.19) 


(4.20) 


The  argument  of  the  MCF,  which  is  the  average  phase  difference  between 
two  received  signals,  is  given  by  the  phase  of  YT  added  to  that  of 


Y^j.  The  characteristics  of  the  individual  coherence  factors  will 


be  analyzed  in  Chapter  5. 


4.3  EXTENSION  TO  SOURCES  SEPARATED  IN  SPACE/TIME 


The  preceding  section  has  derived  the  MCF  for  a  fixed  source 


location.  An  extension  of  the  analysis  to  include  scanning  to  a 


different  location  at  a  later  time  will  introduce  additional  coher¬ 


ence  factors  due  to  the  effects  of  randomness  in  the  scanning  channel. 


The  VLA  system  design  procedure  discussed  in  Sections  1.1  and  3.1 


requires  the  use  of  a  known  beacon  source  upon  which  the  array  can 


initially  focus  due  to  the  unknown  multipath  structure  and  unknown 


phase  of  each  ray  due  to  the  initial  state  of  random  fluctuations. 


The  source-receiver  configuration  for  scanning  is  illustrated 


in  Fig.  4.2.  From  a  beacon  source  at  location  y  and  time  t,  the 


sensor  at  x  receives  a  signal  proportional  to  the  transfer  function 


of  the  channel,  denoted  by 


H(w,x,y,t)  =  A(oj,x,y,t)e^^W,x,y,t:^ 


(4.21) 


and  the  sensor  at  x  +  c,  receives 


H(w,x+S»y,t)  =  A(w,x-K,y,t)ej<Kl0’,H'C’y’t:)  . 


(4.22) 


It  is  desired  to  form  a  VLA  by  focusing  the  receivers  on  the  known 


source  at  y,  t,  and  then  scanning  for  an  unknown  source  at  y-Hi,  t+x, 


Each  receiver  cophasos  for  the  beacon  source  by  using  a  matched 
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1  Random  channel  representation. 


SCAN  LOCATION 
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Fig.  4.2  Source-receiver  configuration  for  scanning 


shift  necessary  for  scanning.  The  quantity 

(w,x,y+n,t+T)  -  <j>’  (w.x+f.y+n.t+T)  (4.30) 

will  be  found  to  be  the  negative  of  the  phase  of  the  coherence  func¬ 
tion. 

The  transfer  functions  for  the  scan  channels  after  cophasing 
for  the  beacon  source  are 

H_<«)  =  H'(w)H*(w)  (4.31) 

m  in  m 

and 

Kn(to)  -  h;(w)H*(o>)  (4.32) 

where  the  subscripts  m  and  n  denote  sensors  at  x  and  respectively 

and  the  prime  denotes  the  scan  channels,  i.e, 

_  Km  -jwT, 

H  (W)  -  H(u,x,y , t)  =  l  A,  e  m  ,  (4.33) 

m  km 

_  Km  -jwT’ 

H '  (w)  =  H(w,x,y+n,t+T)  =  l  A '  e  ,  (4.34) 

m  ktil  km 

_  Kn  -jwT 

H  (w)  -  lI(w,xH»y,t)  =  l  A,  e  Kn  ,  (4.35) 

n  k=l  Kn 

_  Kn  -jwT,* 

H'(w)  =  H(u),x+C,y+n,t+x)  =  l  A '  e  t'n  .  (4.36) 

n  k-1  Kn 


The  KCF  for  scanning  in  space  and  time  is  now 


v  .fn\  _ >>  it  £ _ 


<l«i<“>H*c<o>l2> >s  <|n;(M)nn(«.)  1 2>' 


(A. 37) 


A. 3.1  DISCUSSION  OF  SCANNING  CHANNEL 


The  extension  to  scanning  introduces  the  remaining  type  of 


acoustic  fluctuation,  that  due  to  spatially  varying  multipath  inter¬ 


ference  as  discussed  in  Section  3.2.1. 


In  order  to  determine  this  effect  on  scanning,  the  following 


scan  channel  model  will  be  postulated.  The  scanning  geometry  is 


depicted  in  Fig.  4.3,  in  which  S  is  the  linear  horizontal  scan  dis¬ 


tance  from  the  beacon  to  a  new  source  location.  The  components 


along  the  new  source-receiver  paths  are  designated  x  and  x  ,  and 

m  n 


correspond  to  the  changes  in  source  range  due  to  scanning.  As 


postulated  in  Section  3.1,  S  the  correlation  distance  of 


the  large  scale,  long  period  environmental  fluctuations.  The 


deterministic  multipath  field  in  the  absence  of  the  smaller  scale 


environmental  fluctuations  can  then  be  considered  as  azimuthally 
isotropic  for  a  given  receiver.  As  prescribed  in  Section  3.2.1,  the 


same  rays  are  received  throughout  the  scan  area,  and  the  require¬ 


ment  that  each  ray  describes  a  piano  wave  with  the  same  arrival 


angle  throughout  the  scan  area  is  satisfied  if 


I  W 


P3 
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x  «  R 
m  m 


(A. 38) 


for  each  receiver.  Also,  the  relative  amplitudes  of  the  rays  do 
not  vary  with  changes  in  source  range  due  to  scanning  if  the  above 
condition  is  valid  [11]. 

With  this  realistic  model,  then,  the  total  ray  travel  time  for 
each  channel  can  be  decomposed  as  follows  (the  subscript  m  or  n  is 
implied) : 


beacon  channel 


\ "  \o +  v +  v  ray k 


M.39) 


identified  as 


-  nominal  travel  tir.e  defined  in  Section  4.2. 

-  the  fluctuation  described  in  Section  4.2  which  is  indepen¬ 
dent  and  identically  distributed  among  the  rays  of  the 
beacon  channel,  and  uncorrelated  between  receivers.  It  is 
now  assumed  that  it  is  a  zero-mean  Gaussian  random  process 
with  the  following  characteristics: 


(‘ito'mi.)  "  °-  w  oc  n s 

<0  *l>  r°r  “ji  k> 


(4.40) 


(4.41) 


*  -'iii 


t,£  -  the  correlated  fluctuation  defined  in  Section  4.2. 

scan  channel 


T’k’  -TVo  +  tV  +  t,T  • 

T*k’0  =  Tk'0  +  t’sh'  ’  ray  k' 


(4.42) 
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identified  as 

T*.  -  the  nominal  travel  time  in  scan  channel. 

k'O 

-  the  component  of  T’^iq  which  is  the  nominal  travel  time 
of  ray  k*  in  the  beacon  channel. 

tg^,  -  the  additional  travel  time  in  the  scan  channel  due  to  a 
change,  x,  in  the  source  range,  defined  in  Section  3.2.1 
as 


•4'  ■  ?os6k' 


(4.43) 


'Wk' 


The  0j  ,  were  assumed  to  be  independent  random  samples  from 
the  same  distribution.  An  additional  assumption  is  now 
made  that  the  arrival  angles  are  independent  between  re¬ 
ceivers.  This  is  reasonable,  since  widely  spaced  sensors 
receive  entirely  different  multipath  fields.  The  ray 
arrivals  are  not  plane  waves  across  the  receivers,  and  the 
nominal  travel  times  also  differ  due  to  the  larger  scale 
fluctuations  (note  that  no  restriction  was  made  on  receiver 
spacing  with  respect  to  the  larger  scale  fluctuations; 
due  to  the  large  time  scale,  they  are  frozen  for  all  time 
parameters  of  relevance  in  this  problem,  and  can  therefore 
be  considered  as  deterministic,  contributing  only  to  the 
nominal  travel  times). 

the  fluctuation  described  previously.  However  it  may  now 
be  correlated  with  the  fluctuation  of  ray  k*  in  the  beacon 


*1 
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channel  of  the  same  receiver  if  the  scan  distance  is 
small.  With  the  Gaussian  assumption,  its  characteristics 
can  be  summarized  in  terms  of  rms  values  and  its 
correlation  coefficient  as 


<(tiik’mtMrn')  ‘  °-  kW  °r 

(t'L'm)  ’  4'm-  for  a11  k'-  m' 
^k'nVio)  *  °*  kVi  or  rt"; 
<(tWk'n,tKk'm')  ’  p4»V  tor  a11  k’  •  “• 


(4.44) 

(4.45) 

(4.46) 

(4.47) 


the  fluctuation  which  is  correlated  among  rays  of  the  scan 
channel.  Since  there  may  be  a  correlation  between  receivers, 
this  implies  that  there  may  also  be  a  correlation  between 
the  scan  channel  and  beacon  channel,  since  scan  distance 
will  generally  be  smaller  than  receiver  separation. 


4.3.2  DERIVATION  OF  THE  COHERENCE  FUNCTION 

From  (4.31),  (4.32),  and  (4.37)  the  numerator  of  the  MCF  is 


(h  f/*\  =  /hm*h’*h  \  . 

\  m  ns  \  m  ra  n  n/ 


(4.48) 


Substituting  the  transfer  functions  from  (4.33)  -  (4.36)  yields 
the  expressions 


KA 


H'll  =  £  y  A,  A.  ,  exp— jej (T 1 .  ,  -T,  )  , 

m  m  ,  £  J, ,  km  k'it'i  1  J  k  m  km  ’ 


(4.49) 


LjataML i  fii  :'L--  ^ -  2AJ 
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K  K 
n  n 


B'nHn  ’  l  l.hth'n**  WvnW 


(4.50) 


The  expected  value  in  (4.48)  will  then  be 


(4.51) 


Expanding  the  exponential  into  its  components  gives 
<eSp-jW[Ti,iii-Tta-T-,n+T)ln]V>  = 

<c^!0(t;!k-m-W!><CeXp-3“T  cos6k'm^  ' 
exp  3"(Vn<fTJlnoKexp  <sf‘P  H?  “‘Vn^  ' 

^*p-juCt^-tTm-t^n+tIn)  )>  .  (A. 


i’nO  £n(r 


C4.  v^) 


Denoting  the  three  factors  on  the  above  lines  by  a,  ,,  ,  ao»o  *a 

n  Kin  x»  mn 


CH  H*V)=  (57  A,  A,  ,  a, ,  ,  )(5T  A0  A0I  a,,,  )a  ,  (4,53) 

\  m  n  /  '££,  km  k'm  kk'ra  in  i’n  ii'n  mn  ' 


The  first  component  factor  of  contains  the  phase  due  to  the 
nominal  travel  times  of  channel  ra,  the  second  due  to  uncorrelated 
ray  fluctuations,  and  the  last  due  to  spatial  variations  in  scan¬ 
ning;  the  same  description  applies  to  the  factors  of  a^t£n  for 
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channel  n.  The  factor  Y__  '  contains  the  effect  of  fluctuations 

'Tran 

which  are  correlated  between  channels  m  and  n,  including  the  scan 
channels . 

Consider  now  the  expansion  of  (4.51)  due  to  the  first  factor 
in  (4.52), 


|*(|'Yju^k,mexp"'^tl,^klOra_^!.Jm^^.exp-^tl)^CKklni-Ck'kin^,\cxp-^Wc^COS®k,m^  * 


The  first  expected  value  is 


(4.54) 


<«p-J“<WW>  =  exp  -  ¥2  <C(tW»rW2)> 


exp  -  |jto2D^(S,T)  ,  k’=k 

1  2  (4> * 2  +  $2),  k’#c  , 


Lcxp  -  ^  m 


(4.55) 


in  which  D’ (S,t)  is  the  structure  function  of  t,„  defined  as 
m  ’  Wkrn 

d;(S,t)  =  $1  -  2p(S,T )*  $'  +  i>'2  (4.56) 

id  id  n  m  m 

where  p(S,T)  is  the  correlation  coefficient  from  (4.47).  The 
characteristic  functions  are 


:sJ“5  -(oxp-ju-jl  c°sS 

(4.57) 

k»<u)  '  cxp  -  -i”2,1.2, 

(4.58) 

%  /  \  1  2. ,2 

V,  (w)  =  exp  -  -tm  <>' 

Wra  ‘2m 

(4.59) 

---...JW.  tT’i  I  i  ‘ T  *V‘  -*-2.1 — it 
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so  that  (4.54)  becomes 


cSmtCXP 


“  5“  Dm(S,T)^Akra+  cSrac^Ci^,^k^aaAk,mexp“:Jw(Tk,Om"TkOm)  ‘ 


(4.60) 


An  important  simplification  can  be  made  if  it  can  be  assumed  that 
scan  distance  and  time  are  greater  than  the  correlation  distance 
and  time  of  t^,  i.e.,  S  >  Lq  and  T  >  Tq>  so  that  p(S,t)  K  0.  In 
Chapter  5,  t^k  will  be  identified  with  internal  wave  fluctuations, 
for  which  Lq  «  6.4  km  and  Tq  =  1*6  hr  [4].  Since  the  primary 
interest  of  this  study  is  for  scan  distances  and  times  greater  than 
these  values,  it  will  be  assumed  that  p(S,T)=0.  (This  point  will  be 
discussed  further  in  Chapter  5.)  With  this  simplification  then 

exp  -  |oV(S,T)  =  exp  -  |*>2(*;W)  -  (4*6I) 

so  that  (4.60)  becomes 

CSmCWnCWia^Akm)  ^mO^  =  [cW^\m^HmO^  (4.62) 

which  has  been  factored  into  separate  terms  for  the  beacon  and  scan 

channels  and  where  H  is  the  normalized  transfer  function  of  channel 

raO 

m  in  the  absence  of  fluctuations,  as  defined  in  (4.12). 

The  result  for  the  second  factor  of  (4.52)  is  derived  in  an 
identical  manner.  The  complete  result  for  the  numerator  of  the  MCF  is 
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y  “  f cKm^Akm^  I HmO ^  ^ CSmCWn ^Scm*  I HraO ^  * 

^cWn.^Akn^  ^ HnO ^  ^  ^cSnCwn ^Sen^  ^ HnO ^  ^Tmn  *  (4. 

•4 

The  two  factors  of  the  denominator  of  the  MCF  likewise  have 
identical  derivations.  The  first  factor  is 


(4.64) 


With  the  assumption  made  above  that  p(S,T)  =  0,  the  magnitudes  of 
the  transfer  functions  are  independent  between  the  beacon  channel 
and  scan  channel,  so  that 


<l«jV^^KlV/lHil2\%  •  '  (4.65) 

\  a  /  \  “  /  NU1/ 

The  square  of  the  first  factor  of  (4.65)  is  equation  (4.13), 

<lH/>  -  l4n+<^V2|\o|2-|4»>4  *  «-13> 


and  a  similar  derivation  for  the  scan  channel  yields 


<IH;!2>  -J,4'»+'(|,W2|Hmo|2-^'n,]lCSn,|2c™  -<4-66> 

The  expressions  for  ^  |  | ^  and  <^|h^| are  analogous. 

The  final  result  can  now  be  written  as  a  composite  of  five 
factors, 


v  =  Y  y’Y  Y*  Y 
'Smn  'm'm'n'n  'Tmn 


(4.67) 
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in  which  the  prime  denotes  the  auto-coherence  for  the  scan 
channel.  As  in  (4.14),  the  substitution  K  =  (£a^)  /(£a^)  is 
made  for  each  auto-coherence  factor.  Using  the  envelope  approxi¬ 
mation,  the  results  are 


•m 


K  c2 
m  wm 


1+<V1)cW» 


Wm  'Mm 


(4.68) 


Y’ 

m 


K  |cc  |2c*  2 
m1  Sm'  Wm 


34>, 


Sm 


H 


mO1 


Y '  •  Y* 

TW,Sm  TMm 


J4> 


'Sm 


lCSm'e 


Sm 


(4.69) 


K  c2  n 
n  Wn 


1+(K  -l)c* 
n  Wn 


(4.70) 


K  |c  1 2c * 2 
n 1  Sn 1  Wn 


W'V-Wl'SnlS;2 


“j<i> 


Sn 


CSn  'CSn*G 


Sn 


YTmn  =  ^>exP":}w(tTm"tTm“t:'Tn‘ht:Tn) 


) 


(4.71) 

(4.72) 


The  solution  to  an  extremely  complex  problem  has  been  reduced  to  a 
composite  of  strikingly  simple  factors,  with  no  restrictive 


i 


i 


he 


assumptions  or  approximations.  Equations  (4 - 67)— (4 .72)  are  the 
most  important  results  of  this  work. 

The  first  important  feature  of  this  solution  is  that  it  includes 
the  MCF  without  scanning  developed  in  Section  4.2  as  a  special  case. 
That  solution  is  obtained  by  setting  all  primed  auto-coherences  to 
unity,  and  omitting  the  primed  fluctuations  from  y^^.  (The  result¬ 
ant  phase  of  the  multipath  transfer  functions  does  not  appear  now 
since  the  beacon  is  used  as  a  focus;  also,  the  former  solution 
cannot  be  found  by.  letting  S-K),  since  it  was  assumed  that  S>Lq,  which 
makes  the  scan  channel  and  beacon  channel  independent.) 

The  first  auto-coherence  factor,  equation  (4.68),  is  a  composite 
of  the  effects  of  uncorrelated  ray  fluctuations  and  frequency  select¬ 
ive  fading  in  beacon  channel  m.  Equation  (4.69)  is  the  auto-coherence 
for  scan  channel  m.  The  additional  effect  of  fluctuations  due  to 
spatially  varying  multipath  interference  now  multiplies  the  effect 

of  uncorrelated  fluctuations.  The  phase  of  y*  ,  <j>  ,  is  the  average 

V# 1  om  bin 

phase  difference  between  the  scan  location  and  the  beacon.  It  is  the 
primary  component  of  the  phase  shift  for  receiver  m  which  will  be  re¬ 
quired  for  scanning.  The  auto-coherence  factor  due  to  frequency 
selective  fading  is  the  same  as  that  for  the  beacon  channel,  since  it 
has  been  stipu3ated  that  the  multipath  field  is  azimuthally  isotropic 
over  small  scan  distances.  The  resulting  effect  is  that  the 
extension  to  scanning  has  squared  the  coherence  due  to  frequency 
selective  multipath  interference.  However,  it  will  be  seen  in  Chapter 
5  that  this  has  no  degrading  effect  at  coherent  frequencies. 
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The  auto-coherences  for  channel  n  have  the  same  interpretation 
as  above.  The  last  factor  of  the  MCF  is  the  coherence  due  to 
fluctuations  which  have  some  correlation  among  the  channels,  and 
which  will  be  developed  in  Chapter  5.  The  phase  of  this  term  is  an 
additional  phase  difference  between  channels  m  and  n  required  for 
scanning. 

The  convenient  factorization  of  the  MCF  into  eight  auto-coherence 
functions  and  a  coherence  due  to  correlated  fluctuations  allows  group¬ 
ing  of  terms  to  determine  relative  effects  of  various  combinations. 

To  study  the  relative  contribution  of  scanning  to  coherence,  write 

4m  “  n 


K 


v  Y*  Y 
on  ran'Tmn 


(A. 73) 


and  Y*  can  be  compared  to  y  .  The  relative  contribution  of  each 
'mn  r  ‘mn  . 

receiver  channel  is  similarly  determined  from 


Ysmn  "  <  W  YTmn 


"  Wta  • 


(A. 74) 


by  comparing  Ym  to  Yn»  The  most  important  simplification  is  the 
separation  of  the  effect  of  random  fluctuations  from  that  of  frequency 
selective  multipath  interference  by  writing 
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*Smn  “  (Wi,  i,Wi ,  Sn^Tran  >  ^jUta* 


(\’,sVTM  • 


(4.75) 


The  value  of  this  factorization  is  that,  since  the  effect  of  random¬ 
ness  forms  an  envelope  of  the  MCF  and  is  a  mono tonic ally  decreasing 
function  of  frequency,  it  enables  a  prediction  of  maximum  coherent 
frequencies'  without  knowledge  of  the  particular  multipath  structure 
or  its  frequency  selective  coherence  function,  Y«(w). 


4.4  SUMMARY 

This  chapter  is  the  most  important,  and  the  theory  presented 
provides  the  basis  for  the  rest  of  the  dissertation.  The  theory  of 
the  multipath  coherence  function  has  been  developed  based  upon  a 
formulation  of  the  spectral  coherence  function  in  terms  of  the  random 
multipath  channel  transfer  functions.  This  has  shown  that  the  MCF 
is  independent  of  the  signal  source,  and  depends  only  on  the 
characteristics  of  the  channel.  It  therefore  applies  equally  well 
for  narrow  band  or  broad  band,  random  or  deterministic  signals,  at 
each  frequency  in  the  source  spectrum. 

Due  to  the  stochastic  independence  of  channels,  the  MCF  factors 
conveniently  into  two  auto-coherences.  The  value  of  this  factoriza¬ 
tion  is  that  each  channel  can  be  analyzed  independently,  rather  than 
computing  non-separable  coherences  for  all  pairwise  combinations 


of  receivers. 
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The  MCF  has  been  formulated  to  consider  the  two  types  of 
environmental  fluctuations:  those  which  cause  uncorrelated  ray 
fluctuations  and  those  which  cause  correlated  fluctuations.  The 
MCF  has  been  generalized  to  include  the  latter  type  as  a  cause  of 
acoustic  fluctuations  which  may  be  partially  correlated  between 
receivers. 

The  next  important  development  is. the  envelope  approximation, 
whereby  each  auto-coherence  factors  into  two  coherence  terms,  one 
for  the  effects  of  random  fluctuations  alone,  and  the  other  for 
frequency  selective  multipath  interference.  This  allows  computa¬ 
tion  of  maximum  coherent  frequc.  'ey  independent  of  the  multipath 
configuration. 

The  generalization  of  the  MCF  to  include  the  effects  of  scan¬ 
ning  introduced  another  type  of  acoustic  fluctuation,  that  due  to 
spatially  varying  multipath  interference.  This  fluctuation  was 
accounted  for  by  applying  a  stochastic  model  to  the  ray  arrival 
angles.  Due  to  the  weak  assumption  that  scan  distance  and  time  were 
larger  than  the  corresponding  correlations  of  environmental  fluctua¬ 
tions,  the  MCF  could  again  be  factored  into  separate  coherence 
functions  for  the  scan  channel  and  beacon  channel.  The  resulting 
generalized  MCF  is  a  concise  mathematical  expression  composed  of 
simple  factors  which  allow  any  single  coherence  term  to  be  analyzed 
separately. 

The  remaining  task  to  be  performed  in  Chapter  5  is  the 


specification  of  .the  MCF  parameters  in  terms  of  .real  oceanographic 
fluctuations.  The  parameters  of  environmentally  caused  fluctuations 
will  be  derived  from  the  theory,  of  internal  waves  and  tides,  and  the 
effects  of  both  spatial  and  frequency  selective  multipath  interference 
will  be  determined  from  realistic  models  of  the  underwater  channel. 
However  it  must  be  emphasized  that  the  results  of  this  chapter,  the 
most  important  of  which  are  equations  (4.67)  -  (4.72),  do  not  depend 
upon  the  presently  known  types  of  real  oceanographic  fluctuations  and 
their  actual  stochastic  parameters,  but  only  require  that  they  be 
classified  as  described  in  Section  4.3.1.  Should  future  oceanographic 
developments  provide  an  update  of  the  present  state  of  knowledge,  the 
model  will  still  be  completely  applicable. 
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CHAPTER  5 


■THE  COHERENCE  FUNCTION  IN  TERMS'  OF  THE  OCEANOGRAPHIC  FLUCTUATIONS 

5.1  INTRODUCTION 

In  Chapter  4  the  general  form  of  the  MCF  was  derived  for  beam- 
forming  and  scanning  in  multipath  channels.  The  travel  time 
fluctuations  in'  the  ray  paths  were  defined  in  terms  of  their  general 
stochastic  characteristics,  but  their  parameters  were  not  specified 
in  terms  of  environmental'  fluctuations. 

Chapter  3  identified  the  four  primary  types  of  oceanographic 
fluctuations  which  affect  coherence:  spatial  variations  due  to 
multipath  interference,  internal  waves,  internal  tides,  and  frequency 
selective  multipath  interference.  The  first  three  types  cause 
travel  time  fluctuations  in  the  ray  paths,  and  the  stochastic  para¬ 
meters  of  these  fluctuations  were  specified.  It  now  remains  to 
identify  these  fluctuations  with  those  of  the  MCF  developed  in 
Chapter  4  in  order  to  determine  signal  coherence  in  real  ocean  chan¬ 
nels. 

The  travel  time  of  a  ray  in  the  beacon  channel  was  decomposed  as 

Tk  “  WV  ra*  k  •  w>39> 

In  terms  of  oceanographic  fluctuations  they  are  identified  as 

Tro  ~  nominal  travel  time  affecting  frequency  selective 
multipath  interference. 
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t,„  -  fluctuation  due  to  internal  waves. 

Wk  ' 

-  fluctuation  due  to  internal  tides. 


In  the  scan  channel 


T* 


Tk'0+tWk,+t:T 


^'k’O  Vo^Sk’, 


ray  k*  , 


(A. 42) 


arid  there  is  an  additional  fluctuation, 


tg^,  *r  fluctuation  in  scanning  causing  spatial  variations  due 
to  multipath  interference. 


The  effect  of  each  of  these  fluctuations  on  coherence  will  be  deter¬ 
mined  in  the  following  sections. 


5.2  EFFECT  ON  COHERENCE  OF  OCEANOGRAPHIC  FLUCTUATIONS 

The  system  geometry  for  scanning  away  from  a  beacon  using  a 
two-receiver  array  was  described  in  Section  4.3.1  and  illustrated  in 
Fig,  4.3.  The  purpose  of  this  section  is  to  determine  the  MCF 


mn 


YMnYNraYW, SmYMmYtfnYNnYW, SnYMnYTmn 


(5.1) 


where  the  individual  auto-coherence  factors  were  defined  in  (4.68)- 
(4.72).  In  terms  of  oceanographic  fluctuations  they  are  now 
identified  as 

Y-.  -  effect  of  internal  waves  in  channel  from  beacon  to 
Wra 

receiver  m. 

YMm  ”  cY^cct  °f  frequency  selective  multipath  interference  in 


Y* 

'tf,Sm 

YTmn 


beacon  channel  to  receiver  m. 

-  effect  of  internal  waves  and  spatial  variations  due 

to  multipath  interference  in  scan  channel  to  receiver  m. 

-  effect  of  internal  tides  in  beacon  channels  and  scan 
channels  to  both  receivers  m  and  n. 


The  remaining  factors  in  (5.1)  have  corresponding  definitions  for 
receiver  n  or  for  the  scan  channel  (denoted  by  a  prime) .  The  con¬ 
tribution  of  each  type  of  fluctuation  to  the  MCF  and  its  relative 
importance  will  now  be  determined  in  terms  of  its  respective  auto¬ 
coherence  factor. 


5.2.1  INTERNAL 'WAVES 


A  basic  premise  of-  this  work  has  been  that  the  receivers  are 
separated  by  such  large  distances  that  travel  time  fluctuations 
induced  by  internal  waves  are  independent  between  them.  In  Section 
4.3.2  it  was  further  assumed  that  horizontal  scan  distance,  S,  and 
scan  time,  T,  are  larger-  than  the  corresponding  correlation  distance 
and  time  of  the  fluctuations,  so  that  the  fluctuations  in  the  scan 
channel  are  independent  of  those  in  the  beacon  channel.  In  Section 
3; 2. 2  the  correlation  coefficient  was  given  as 


P(S,t) 


(3.12) 


and  is  illustrated  in  Fig.  5.1.  From  this  equation  the  scan  dis¬ 
tance  S  for  which  the  beacon  and  scan  channels  are  independent  can 


be  determined  for  a  given  time  T  from,  initial  focus  on  the  beacon 
The  auto-coherence  due  to  internal  waves  in  each  beacon 


.channel  is  of  the  form 


VT 


Kc 


W 


l+(K~l)c 


W 


(5.2) 


In  Section  4.3.2  the  characteristic  function  was  shown  to  be 


cw  “  exp(-  |  wV)  (5.3) 

as  is  shown  in  Fig.  5.2  as  a  function  of  f$.  The  mean  square  travel 
time  fluctuations  were  given  as 

4>2  **  (3.4x10  ®sec2km  X)R,  steep  ray;  (3.9) 

4>2  *=  (6,8x10  ®sec2km  ^)R,  axis  ray,  (3.10) 

and  are  shown  in  Fig.  5.3. 

With  these  equations  the  auto-coherence  due  to  internal  waves 
for  each  channel  can  be  computed  as  a  function  of  acoustic  frequency 
and  the  range  to  the  beacon  from  each  receiver.  Fig.  5.4  illustrates 
a  typical  variation  of  y^  with  beacon  range,  and  the  attenuation  with 
acoustic  frequency  is  depicted  in  Fig.  5.5.  Both  computations  assume' 
steep  rays  using  (3.10)  and  the  ray  parameter  is  K  =  4. 


5.2.2  SPATIAL  VARIATIONS  DUE  TO  MULTIPATH  INTERFERENCE 

The  effect  of  spatial  variations  due  to  scanning  for  each  channel 


SCAN  UME,  T  (HRS) 


TRAVEL  TIKE  FLUCTUATION,  $  (MSEC) 
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is  determined  from 


Y*  «s 

Tw,s 


*k|2c’v 


l+(K-l)|cs|2c'f2 


(5.4) 


The  characteristic  function  for  the  spatial  variations, 


cs  "  lcsle 


(5.5) 


was  developed  in  Section  3.2.1  and  can  be  written  in  terms  of  the 


wavenumber,  kg,  as 


r  2*1  f  *" 

ll+(k0>:a2)  exp  j  -kQx+  -|tan  1(kQxa2)  . 


(5.6) 


The  magnitude  of  cc  consists  of  the  first  factor.  Fig.  5.6  shows  the 


variation  of  |c_|  with  |x|/X  for  characteristic  values  of  the  ray 


spread,  0. 


In  Section  5.2.3  it  is  shown  that  internal  tides  have  no  effect 


on  average  signal  phase.  Therefore  the  term  4>  is  the  total  average 


phase  change  for  one  receiver  channel  due  to  scanning  away  from  the 


beacon.  In  (5.6)  it  is  seen  to  consist  of  two  terms.  The  first 


term  is  the  linear  component,  -kgX.  The  second  component  is  due  to 


the  ray  spread.  Note  that  <*><.(-x)  =  -!{>q (x)..  The  phase  with  the 
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Fig.  5.5  Frequency  variation  of  internal  wave  auto-coherence. 
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Fig.  5.6  Characteristic  function  for  spatial  multipath  interfer- 


linear  component  removed  is  illustrated  as  a  function  of  x/X  in 
Fig.  5.7  for  characteristic  ray  spreads. 


The  coherence  y'  _  also  includes  the  effect  of  internal  waves 

W  f  o 


in  the  scan  channel.  For  the  purpose  of  comparison  with  Yw>  the 


characteristic  function  c^  is  set  equal  to  unity  and 


Iyc 


KjcJ‘ 


1+CK-Dlcgl4 


(5.7) 


is  computed.  Assuming  a  ray  spread  a  B  2°,  the  variation  of  | Yg I 
is  illustrated  in  Fig.  5.8  as  a  function  of  |x|,  and  in  Fig.  5.9 
as  a  function  of  frequency.  Note  the  larger  rate  of  attenuation 
of  | Ysl  with  range  and  frequency  compared  to  that  of  y^  in  Figs.  5. A 
and  5.5.  This  indicates  that  for  a  given  increase  in  range  due 


to  scanning,  the  decrease  in  | J  is  much  more  severe,  than  the 
corresponding  decrease  in  y^,  and  is  the  limiting  factor  in  scanning 
ability.  Since  x  is  the  change  in  range  due  to  scanning,  it  can 
also  be  concluded  that  the  maximum  limitation  on  scanning  is  in  the 
direction  of  the  propagation  path  from  beacon  to  receiver.  In  a 
direction  perpendicular  to  this  path  the  change  in  range  is  much 
less  so  that  there  is  less  limitation  on  scanning. 
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Fig.  5.7  Average  phase  variation  due  to  spatial  multipath 
interference. 
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Fig.  5.8  Range  variation  of  auto-cohcrcnce  due  to  spatial 
multipath  interference. 
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5.2.3  INTERNAL  TIDES 

The  effect  of  internal  tides  on  coherence  is  described  by 

YIbn  "  <cxp-3“<til,rtTm‘tin+';l„)> 

*  ^xp-jwAt^  .  (5.8) 


The  travel  time  of  an  axis  ray  from  source  to  receiver  in  the 
presence  of  internal  tides  was  derived  in  Section  3.2.3  as 


1  -  2Ac0 


sin  (w_t-k_Rcos'{>  /2  )- 


sin(kTRcos({)/2) 

(kTRcos;j>/2) 


(3.18) 


where  =  R/cq.  The  travel  time  fluctuation  due  to  internal  tides 
is  the  same  for  each  ray  [1]  so  that  the  results  for  an  axis  ray  are 
used. 


Consider  the  simplified  source-receiver  configuration  shown  in 
Fig.  5.10.  Two  sensors  separated  by  a  distance  Rc  are  located  on  a 
baseline  perpendicular  to  the  direction  of  internal  tide  propagation 
(e.g.,  on  a  continental  shelf).  A  beacon  is  located  equidistant  from 
* two  receivers  at  a  range  Rq.  At  time  t,  the  travel  times  to  the 
two  sensors  are 
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(5.9) 


0 


30 


60  .90 

FREQUENCY,  f  (HZ) 


120 


Fig.  5.9  Frequency  variation  of  auto-coherence  due  to  spatial 
multipath  interference. 
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—  sin(wTt-kTRncos<J)M/2) 
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(5.10) 


The  travel  times  from  the  scan  location  at  (S,0),  at  some  later  time 
t  '+  T  are 
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(5.11) 


.  and 
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sin[uT(t+T)-kTR^cos^/2] 


sin (k  R'cosfJ)* /2) 
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(5.12) 


The  travel  time  fluctuations  due  to  the  internal  tide  are 
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But  tTra  *  tTn*  so  that 


(5.14) 


It  is  also  true  that 


R’cosi'  =  R'costj)'  .=  R.+SsinO 
n  n  m  in  0 


so 


(t+T)-kT(R0+SsinO)] 


sin[kT(P0+Ssin6)/2] 
^Ik^(Ro+Ssin0)/2]  * 

(5.15) 


Here  the  quantity  AR  %t  is  the  range  difference 

mu 


AR  «  R'-R' 
mn  m  n 


(5.16) 
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the  effect  on  coherence  is  given  by  the  factor 


(5.17) 


which  will  denote  an  ensemble  average  over  all  time  of  initial  focus 
on  the  beacon  0  ^  t  <■•  2tt/w  ,  i.e. 


(5.18) 
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Writing  the  phase  as 


wAtj  =  asinfw^Ct-J-T)-^ 


(5.19) 


then 


-jwAt 


l>  -  ^ 


-jctsin[toT(t+r)-^] 


(5.20) 


JQ(a) 


in  which  Jq  is  the  zero  order  Bessel  function.  The  complete  effect 
on  coherence  between  channels  tn  and  n,  due  to  internal  tides,  is 
therefore  given  by  the  expression 


’2Ac0 \  [ AR J  sin tlcT (RQ+Ssine ) /2 ] 
C0  j(  c0  }  [kT(RotSsin0)/2] 


(5.21) 


The  only  assumption  which  has  been  made  in  this  derivation  is  that 
Acq/cq  «  1  (a  characteristic  value  for  Ac^/c^  due  to  internal  tides 

-5 

is  10  ).  Since  is  real,  it  makes  no  contribution  to  the 

average  phase  difference  between  the  signals. 

It  is  important  to  analyze  the  physical  significance  of  this  re 
suit.  For  this  purpose,  assume  that  the  scan  distance  S  <<  R^.  It 
can  then  be  shown  that 


AR  ~ 


2Scos0 


-  (5.22) 


The  coherence  then  becomes 


Y  **  J 
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sin(kTRQ/2) 
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(5.23) 


First,  there  is' a  noticeable  absence  of  dependence  on  the  time 
difference,  r.  This  is  explained  by  the  fact  that  the  bulk  time  de¬ 
lays  are  equal  for  the  first  source  location.  If  they  were  not 
chosen  to  be  equal,  the  mathematics  would  become  unwieldy,  but  it 
can  be  shown  that,  in  general,  the  effect  on  coherence  would  be  a 
dependence  cn  a  sinusoidal  function  of  uyr.  This  would  cause  the 
coherence  to  oscillate  between  unity  and  some  minimum  value  deter¬ 
mined'  by  tlie  other  parameters.  The  configuration  considered  here 
corresponds  to  the  minimum  vajue. 

The  manner  of  dependence  of  Y_  on  the  quantities  to  and 

Tran 

(&cQ/cQ)  is  obvious.  The  effect  of  the  quantity  ScosO  is  interest¬ 
ing:  the  coherence  depends  primarily  on  the  component  of  scan  dis¬ 
tance  perpendicular  to  the  tide  normal.  This  is  consistent  with  the 
previous  observation  that  the  maximum  effect  on  phase  fluctuations 
is  when  ray  propagation  is  perpendicular  to  the  direction  of  the 


K 


As  the  quantity  R-/R-  becomes  large,  the  difference  in  travel 
U  o 


time  between  the  two  sensor  channels  for  a;  constant  scan  distance 


S  becomes  small,  causing  coherence  to  increase.  Likewise,  as 


Rq/X^,  increases,  coherence  increases.  The  explanation  for  this  is 


the  fact  that,  since  R^  is  the  component  of  the  ray  paths  in  the 


direction  of  the  internal  tide  propagation,  as  R^/X^  becomes  large 


the  ray  travels  through  a  larger  number  of  periods  of  the  internal 


tide,  and  the  positive  and  negative  variations  of  the  sound  speed 


variations  tend,  to  average  out  to  zero.  Note  that,  when  the  ray 


has  travelled  through  an  integer  number  of  periods  of  the  internal 


tide,  the  sound  speed  variations- are  completely  cancelled  out,  and 


coherence  becomes  unity,  i.e.. 


sin(k  R  /2)  R 

(Von) 


(5.24) 


This  is,  of  course,  exactly  true  only  for  axial  rays  as  considered 


here;  however  it  can  be  concluded  that,  in  general,  coherence  is 


greater  when  acoustic  propagation  is  in  the  direction  of  the  internal 


tide. 


The  coherence,  Y^nin»  is  plotted  in  Fig.  5.11  r.s  a  function  of- 


w6T  where 


(2Ac AR  \  sin(k  *(Rrt+Ssin9)/2]  ; 
_ 0  I  _ mn  i  _  i  0 _ _ _ _ 

cQ  1 1  cQ  I  [kx(R0+SsinO)/2] 


(5.25) 


is  the  travel  time  variation  due  to  the  internal  tide.  In  the  deep 


i  liVit 
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ocean  it  has  been  found  that  the  4m  internal  tide  is  predominant. 

For  a  typical  sound  speed  profile  with  the  sound  axis  at  a  depth 
of  1200m,  reference  [2]  gives  the  sound  speed  variation  as 
ACq  *■  .06m/sec  for  *  1489.55  m/sec.,  so  that  Acq/cq  *»  4.03  x  10~^. 
Fig.  5.12  shews  the  corresponding  variation  of  6T  with  scan  distance 
S,  for  0  «  0®  and  Rg  *  150  km,  and  for  selected  values  of  range  Rq. 
For  other  amplitudes  of  the  internal  tide,  the  appropriate  value  of 
ACq/cq  should  be  substituted  in  (5.25). 

Based  on  the  results  derived  here,  it  will  be  shown  in  Section 
5.3  that  internal  tides  have  a  negligible  effect  on  coherence  com¬ 
pared  to  internal  waves  and  spatial  multipath  interference. 

5.2.4  ■  FREQUENCY  SELECTIVE  MULTIPATH  INTERFERENCE 

The  effect  of  frequency  selective  muLtipath  interference  on  the 
MCF  is  given  by 


'M 


YMmYMraYMnYMn 


(5.26) 


The  individual  auto-coherence  has  the  form 


YMm 


(5.27) 


where  II-  (w)  is  the  normalized  transfer  function  of  the  channel  in 
Urn 

the  absence  .of  random  fluctuations.  The  effect  of  Yw  (w)  on  the  total 

Mo. 

coherence  is  best  determined  by  assuming  K  rays  which  arrive  with 
equal  time  spacings  and  equal  amplitudes.  Following  [3],  the  rays 
arrive  ov  r  an  interval  of  time  T,  which  is  the  time  spread  of  the 

O 
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channel*  end  Tq  is  the  bulk  time  delay  of  the  channel  as  depicted 
in  Fig.  5.13.  The  auto-coherence  is  easily  found  to  be 

sin(wT  /2) 

YMra^  “  Ksin(u)Tg/2K)  •  (5.28) 

The  periodic  lobe  structure  of  this  function  determines  the  actual 
coherent  frequencies,  i.e.,  the  frequencies  of  the  primary  maxima 
of  the  structure  where 

w  T  “  2ivsrK  , 
n  S  ’ 

fn  *=  ~  ,  n  »  1,  2 .  (5.29) 

S 

As  the  time  spread  of  the  channel  increases  for  a  given  number  of 

rays,  there  is  an  increasing  number  of  coherent  frequencies  in  a 

given  bandwidth.  This  is  the  case  for  ■isv.creasing  source  range. 

Also  for  a  constant  value  of  Tg,  the  spacing  between  coherent 

frequencies  increases  as  K  increases,  as  would  occur  upon  entering 

a  convergence  zone.  Also,  if  the  time  spread  is  proportional  to  the 

number  of  rays,  the  location  of  coherent  frequencies  does  not  change. 

The  coherence  bandwidth  centered  on  f  is  determined  by 

n 

Af  =  ~  .  (5.30) 

S 

For  long  range  propagation,  Tg  is  on  the  order  of  seconds,  so  that 

\ 


Af  is.  generally  less  than  1  Hz.  Although  exact  only  for  ray  arrivals 
which  are  equally  spaced  on  the  time  axis,  (5.29)  and  (5.30)  are 
reasonable  order  of  magnitude  estimates  for  arbitrary  multipath 


fields,,  given  K  and  Tg(see  footnote  to  (4.1)). 

The  effect  of  scanning  is  to  square  the  auto-coherence  factor 
for  each  channel  so  that 


^Mm^Mrn 


Ho»c“> 


(5.31) 


The  effect  of  squaring  this  factor  is  to  narrow  the  peaks  and  widen 
the  nulls  of  the  interference  pattern  causing  an  effective  decrease 
in  the  coherent  bandwidth  to 

Af  -  ~~  .  (5.32) 

However  there  is  no  effect  exactly  at  the  peaks  of  the  pattern,  and 
coherent  frequencies  will  remain  the  same, 

. The  total  coherence  is 

Y„-  Il'0n(“)|2|,>0n<“)|2  '  <5-33> 

Since  receivers  spaced  by  large  distances  may  receive  entirely  dif¬ 
ferent  multipath  fields,  the  resulting  coherent  or  partially  coherent 
frequencies  must  be  computed  by  multiplying  the  coherence  factors  for 
each  receiver  as  indicated  by.  (5.33).  However  if  the  sdnsors  receive 
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Identical  multipath  fields,  the  coherent  frequencies  remain  the 
same,  but  the  coherent  bandwidth  is  reduced  to 

Af  -  4-  •  (5.34) 

S 

An  example  of  the  variation  of  y^  with  frequency  for  identical 
multipath  fields  with  K  =  4  and  Tg  =  4  sec  is  illustrated  in 
Fig.  5.14.  For  these  parameters  it  is  found  that  f  =  1  Hz, 

2  Hz,  ...,  and  Af  =  .0625  Hz. 

5.3  THE  C0J5PLETE  MULTIPATH  COHERENCE  FUNCTION 

The  previous  sections  have  presented  the  effects  of  the  individual 
oceanographic  fluctuations  on  their  respective  auto-coherence  factors. 
It  now  remains  to  compare  the  various  effects  and  to  determine  their 
combined  effect  on  coherence.  A  summary  is  then  given  with  respect 
to  the  application  of  these  results  to  the  computation  of  coherence. 


5.3.1  COMBINED  EFFECTS  ON  COHERENCE 

In  Section  4.3.2  the  MCF  was  factored  into  an  envelope  due  to 
travel  time  fluctuations  and  a  coherence  terra  due  to  frequency  select¬ 
ive  multipath  interference  which  was  written  as 


Smn 


(YW,SY'PYM 


(4.75) 


Since  the  first  factors,  due  to  random  travel  time  fluctuations 
alone,  decrease  monotonically  with  frequency,  it  is  appropriate  that 


i 

l 
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they  be  considered  separately.  The  factor  y^  containing  the 

frequency  selective  effects  on  coherence  gives  the  coherent 

.frequencies  for  which  y^  =  1. 

The  four  components  of  the  factor  y  „  all  have  the  same 

W  >  o 

functional  form  written  as 


Y(c) 


l+(K-l)c2 


(5.35) 


and  is  shown  in  Fig.  5.15  for  various  values  of  K. 

There  is  a  subtle  dependence  on  the  ray  parameter  K  (equal  to 
the  number  of  rays  when  they  have  equal  amplitudes).  Since  this 
form  was  obtained  from  the  envelope  approximation  in  Section  4.3.2, 
each  corresponding  auto-coherence  has  a  companion  factor  due  to 
frequency  selective  multipath  interference.  Consider  a  coherent 
frequency  of  this  factor  obtained  from  the  equal  time  spacing  formula¬ 
tion,  say  f^,  and  keep  it  constant  while  increasing  K  so  that  the 
corresponding  factor  of  y^  equals  unity.  Since  K  satisfies 


K 


f  T 
n  S 

n 


(5.36) 


this  can  be  accomplished  by  allowing  K  to  increase  by  increasing  Tg. 
From  (5.35)  it  can  be  seen  that  the  auto-coherence  then  increases 
as  K  increases.  The  explanation  for  this  is  that  coherence  is 
.primarily  determined  by  the  variations  of  the  resultant  phase  of  a 
single  frequency  component  of  each  multipath  signal.  For  a  given 
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phase  variation  in  the  individual  rays*  as  the  number  of;  independent 
rays  increases,  the  variation  in  the  resultant  phase  decreases.  This 
phenomenon  has  actually  been  observed  in  convergence  zones,  i.e.» 
where  many  ray  paths  converge  in  a  focal  zone  [4J . 

The  above  effect  must  be  carefully  considered  in  analyzing  the 
effects  of  spatial  variations  due  to  multipath  interference.  Although 
Y*  c  will  increase  with  K  when  the  other  parameters  are  held  constant, 
the  ray  spread  O  may  also  increase  due  to  the  increase  in  the  number 
of  rays,  and  this  will  'cause  a  decrease  in  y'  q.  The  relationship 
between  K  and  should  therefore  be  considered  in  computations  ot 

YW,S* 

For  the  purpose  of  comparing  the  various  effects  on  coherence, 
the  simplified  geometry'  of  Fig.  5.1Q  will  be  used.  Each  auto-coherence 
term  in  (5.L)  due  to  random  travel  time  fluctuations  was  computed  as 
a  function  of  scan  distance  for  9=0  (perpendicular  to  Rq)  and 
0  »  ±90°  (parallel  to  Rq),  which  are  the  approximate  directions  of  ex¬ 
trema  of  the  variations  due  to  internal  tides  and  spatial  multipath 
interference.  The  ranges  used  are  Rq  =  250  km  and  Rg  =  150  km,  the 
multipath  parameters  are  K  =  16 .and  a  =  2°,  and  the  acoustic-  - 
frequency  is  f  =  50  }Iz.  Figs.  5.16-5.18  show  the  results  for  scan 
distances  up  to  50  km. 

Fig.  5.17  illustrates  the  results  for  0=0,  which  is  the  direc- 
|  tion  of  the  maximum  effect  of  internal  tides,  and  the  approximate 

j  minimum  of  spatial  fluctuations.  The  solid  lines  are  the  approximate 

f 

region  of  validity  of  the  assumption  of  independence  between  scan 

j. 

L 
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channel  and  beacon  channel.  The  dashed  lines  are  extrapolated  to 

give  the  proper  coherence  of  unity  at  S  =  0. 

The  highest  coherence  factor  is  which  remains  at  unity 

throughout  the  entire  scan  distance.  It  was  shown  that  the  effect  of 

internal  tides  decreases  with  increasing  range  Rq,  while  all  other 

effects  increase.  The  conclusion  is  that  internal  tides  have  a 

negligible  effect  on  coherence  for  long  range  propagation  and  for 

scan  distances  of  this  magnitude,  and  henceforth  they  may  be.  ignored. 

This  result  removes  any  restrictions  on  the  system  configuration  or 

its  orientation  with  respect  to  the  direction  of  internal  tide 

» 

propagation  as  in  Section  5.2.3.  Furthermore  it  was  shown  in  Section 
5.2.3  that  internal  tides  have  no  effect  on  coherence  phase. 

Next  in  value  are  the  auto-coherences  due  to  internal  waves  in 
the  beacon  channels,  which  are  equal  due  to  system  geometry  and  do 
not  vary  with  scan  distance. 

The  auto-coherences  Y*  „  and  Y,\  _  due  to  the  combined  effect 

W,  Sm  1 W,  Sn 

of  internal  waves  and  spatial  multipath  interference  in  the  scan  chan¬ 
nels  have  the  largest  effect  on  coherence.  However  0  *=  0°  is  the 
direction  of  the  approximate  minimum  effect  of  the  spatial  variations, 
due  to  smaller  changes  in  range,  so  that  total  coherence  should  be 

higher  in  this  direction.  The  difference  in  the  values  of  y*  „ 

w,  Sm 

and  y^,  gn  is  due  to  differences  in  scanning  ranges.  The  composite 
MCF,  Yginn  is  largest  in  the  direction  9  ~  0°  so  that  this  is  the 
direction  of  largest  scan  distance  for  a  constant  coherence. 


Fig.  5.16  demonstrates  coherence  for  0  =  +90°.  The  coherence 
factors  y_  ,  yt.  and  Yt.  are  the  same  as  in  Fig.  5.17.  The  increase 
in  range  for  a  given  scan  distance  is  the  greatest  in  this  direction. 
The  effects  of  both  the  spatial  variations  and  internal  waves  there¬ 
fore  are  greater  than  in  any  other  direction  and  the  auto-coherences 
Sm  and  **W  Sn  (e<*ua*  by  symmetry)  attain  their  absolute  minimum 

values.  The  MCF  y.,  is  minimum  in  the  direction  0  =  +90°  and  scan- 
'Smn 

nirig  ability  is  consequently  the  most  limited. 

Fig.  5..18  shows  the  effect  of  scanning  in.  the  direction  0  =  -90°. 
The  effect  of  spatial  variations  is  approximately  the  same  as 
0  **  +90°  for  a  given  S,  but  since  range  from  the  receivers  to  the 
scan  location  is  decreasing,  the  effect  of  internal  waves  is  somewhat 
less  than  G  =  +90°.  This  accounts  for  the  slightly  higher  values  of 
*^W  Sm  an(*  Sn  <:aus*nS  a  slight  increase  in  the  MCF,  Yg^*  However 
for  scan  distances  of  the  magnitude  considered  here,  the  difference 
in  the  MCF  between  0  =  +90°  and  0  =  -90°  is  raiuimal  and  scanning 
ability  is  approximately  the  same  in  these  directions. 

The  average  signal  phase  for  each  receiver  channel  varies  as  a 
function  of  scan  distance  according  to  the  change  in  source  range. 
Negative  values  of  phase  correspond  to  increases  in  source-receiver 
range  relative  to  the  beacon,  and  positive  values  indicate 
decreases  in  range.  The  primary  component  of  the  phase  is  the  linear 
variation  h^x.  It  can  be  seen  from  (5,6)  that  for  large  values  of 
k^Jxlo^,  the  magnitude  of  the  phase  is.  approximately  |<j>g|  =  k^lxj  - 
The  decrease  in  Rey^^  when  scanning  with  the  plane  wave  phase  k^x 
rather  than  <{>t,  is  ns  large  as  l-cus^r  -  .293. 


SCAN  DISTANCE,  S 

Fig.  5.16  Comparison  oC,nuto~cohorencc 
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5.3.2  COMPUTATION ' OF  THE  COHERENCE-  FUNCTION 

The.  purpose  of  this  section  is  to  summarize  the  procedure  for 
computation  of  the  MOF.  -It  has  been  emphasized  that  the  MCF  can  be 
computed  for  all  receiver-  pairs  by  the  computation  of  only  the  auto¬ 
coherence  for  each  receiver.  The  following  outline  gives  the  pro¬ 
cedure  for  computation  of  receiver  auto-coherence,  and  the  MCF  for 
each  receiver  pair  is- computed  by  multiplying  their  auto-coherences. 
Procedure 

1,  For  a  given  sound  speed  profile,  beacon  depth,  receiver 
depth,  range  R  to  beacon,  and  frequency  f,  compute  the 
number  of  ray  arrivals  K' ,  relative  pressure  amplitudes, 
A^,  travel  times,  TkQ,  and  arrival  angles,  0^  (usually 
from  a  ray  tracing  program). 

2. 


3. 

A. 


Compute  the  ray  parameter 


and  estimate  the  rms  ray  arrival  angle  from 


(5.37) 


“  K*  „ 

«  b  l  K 


K1  L  V  L- 

K  k~l  k 


(5.38) 


For  given  scan  location  determine  new  range  to  receiver, 
R* ,  and  compute  x  =  R’-R. 

Determine  <5>  =  <!>(R)  and  <5>'  =  ^(R* )  from  Fig.  5.3  or  from 


(3.9)  and  (3.10)  for  characteristic  ray  type  in  channel. 


124 


•  Find  c„  =  ctf(f4>)  and  c '  =  cT,(f$’)  -from  Fig.  5.2  or 
from  (5.3). 

5.  Determine  |cg|  =• |cg(x/X,C) j  from  Fig.  5.6 
or  from  (5.6). 

6.  From  Fig.  5.15  or  from  (5.35)  compute  y(c^)  and  y(c^|cg|) 
for  the  value  of  K  found  in  (5.37). 

7.  Determine  the  phase  <*>g  =  $g(x/X,c)  from  Fig.  5.7  or 
from.  (5.6). 

8.  If  f  is  a  coherent  frequency  (y^  =  1),  the  complete 
auto-coherence  is 

’  .  i4»s 

Y(cw)Y(cw|cg|)e 

,  Coherent  frequencies  are  determined  from  .  ‘ 

K’  -jwT  /  K’ 

1  (oj)  |  =  |  I  \e  k0|  /  l  ^  =1.  ,  .  (5.39) 

k—1  /  k-1 

Thia  above  procedure  is  performed  for  each  receiver  -channel.  In 
terms  of  these  auto-coherences  for  N  receivers,  y  ,  ra  =  1,2,  ....  N, 
the  MCF  for  each  pair  of  receivers  is 

/  j1 

YSmn  =  YmYn’  ra*  n  =  1>2»  ••••  N,  m  ^  n.  (5.40) 

* 

5.4  SUMMARY  .  . 

This  chapter  has  presented  the  MCF  in  terms  of  real  oceanographic 
fluctuations,  and  has  compared  the  effect  of  each  type  of  fluctuation. 
The  condition  for  which  the  scan  channels  arc  independent  for 


125 


internal  wave  fluctuations  was  shown  to  depend  upon  scan  distance, 

S,  and  scan  time,  T.  For  T  *=  0,  p(S,t)  *>  ,5  for  S  »  6.4  km,  so  by 
restricting  the  analysis  to  S  >  6.4  km  the  channels  can  be  considered 
to  be  independent,  and  the  time  dependence  can  also  be  ignored.  The 
coherence  due  to  internal  wave  fluctuations  was  shown  to  decrease 
with  both  range  and  frequency  and  to  increase  with  the  number  of  rays. 

Spatial  fluctuations  due  to  multipath  interference  were  shown 
to  have  the  most  severe  effect  on  scanning,  and  their  effect  is  com¬ 
bined  with  that  of  internal  waves  in  the  scan  channel.  Their  effect 
on  coherence  depends  upon  a  difference  in  range  to  the  receiver  be¬ 
tween  the  beacon  and  the  scan  location.  This  implies  that  the  maximum 
scanning  ability  is  generally  perpendicular  to  the  direction  from 
receiver  to  beacon.  Scanning  is  much  more  limited  in  the  parallel 
direction.  The  coherence  decreases  with  increasing  angular  ray 
spread,  frequency,  and  scan  distance;  it  increases  with  an  increasing 
number  of  rays  within  the  same  spread  of  arrival  angles.  The  total 
average  signal  phase  to  the  scan  location  is  determined  by  the  spatial 
fluctuations  and  each,  receiver  uses  this  as  an  average  phase  shift 
for  scanning. 

The  coherence  due  to  internal  tides  decreases  with  increasing 
scan  distance  and  frequency,  but  increases  with  range.  However,  the 
effect  of  internal  tides  is  negligible  compared  to  the  other  effects 
for  the  scan  distances,  ranges,  and  frequencies  of  interest  here. 
Internal  tides  also  nave  no  effect  on  average  signal  phase. 


7"™  -V,; ' K  «< 


126 


i  ; 


The  above  effects  form  a  roonotonically  decreasing  coherence 
envelope  of  the  effects  of  frequency  selective  multipath  interference. 
This  latter  effect  depends  upon  the  constructive  and  destructive  inter¬ 
ference  of  the  rays  as  frequency  varies.  It  can  be  stated  in  general 
that  the  spacing  between  coherent  frequencies  decreases  with  increas¬ 
ing  time  spread  and  decreasing  number  of  rays,  and  that  the  coherence 
bandwidth  (about  a  coherent  frequency)  decreases  with  increasing  time 
spread.  However,  the  exact  interference  pattern  must  be  computed 
from  the  ray  amplitudes  and  travel  times.  The  dependence  of  the 

auto-coherence  on  Yw  is  determined  by  the  location  of  the  coherent 

M 

frequencies.  Rather  than  compute  for  an  arbitrary  frequency  (since 
Y^j  may  be  low  due  to  destructive  interference),  the  approach  taken  has 
been  tc  assume  location  at  a  coherent  frequency  so  that  Yjjb  1.  Since 
there  generally  will  be  small  spacings  between  coherent  frequencies, 
the  preferred  approach  is  to  determine  coherent  frequencies  from  the 
exact  multipath  summation,  and  to  assume  that  the  signal  bandwidth  is 
large  enough  to  include  at  least  one  coherent  frequency.  This  fre¬ 
quency  is  then  used  for  computation  of  the  coherence  envelope.  This 
subject  will  be  discussed  further  in  Chapter  6, 

The  complete  auto-coherence  can  be  computed  simply  from  the  equa¬ 
tions  and  figures  given  in  this  chapter.  With  the  aid  of  a  ray  tracing 
computer  program  or  other  data,  the  procedure  of  Section  5.3.2  can 
be  used  to  rredict  the  MCF. 
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CHAPTER  6 


APPLICATION  TO  A  SUPERARRAY  SYSTEM  DESIGN 

6.1  SYSTEM  DESIGN  APPROACH 

Array  processing  was  discussed  in  Chapter  2  and  the  VLA  was  com¬ 
pared  with  a  conventional  array.  In  particular,  a  VLA  of  conventional 
subarrays  was  discussed,  and  its  advantages  with  respect  to  gain  and 
beam  pattern  were  emphasized.  In  Chapter  5  the  final  formulation  of 
the  MCF  was  presented  in  terms  of  known  oceanographic  fluctuations. 

The  purpose  of  this  chapter  is  to  apply  the  results  derived  from  the 
MCF  to  a  VLA  of  subarrays. 

Consider  a  system  of  Ny  widely  spaced  conventional  subarrays, 
each  of  whicti  has  Nc  sensors  and  bearovidth  A0,,.  A  beacon  is  placed  at 
B  at  the  range  Rq,  and  the  beam  of  each  subarray  is  scanned  to  this  loca¬ 
tion  as  shown  in  Fig.  6.1(a).  The  beacon  radiates  a  waveform  which 
enables  each  receiver  to  measure  the  impulse  response  of  the  channel. 
With  this  information  each  receiver  then  focuses  on  the  beacon  as  des¬ 
cribed  in  Section  4.3.  The  pattern  of  the  system  then  changes  from  the 
independent  patterns  of  the  subarrays  to  the  near  field  pattern  of  a 
VLA  with  a  high  resolution,  coherent  focus  on  the  beacon  as  shown  in 
Fig,  6.1(b).  Any  ambiguities  in  the  VLA  pattern’ are  limited  to  the 
original  area  of  intersection  of  the  subarray  beams. 

At  the  beacon  the  coherence  is  unity  for  all  subarray  pairs,  so 


129 


that  the  VLA  gain  attains  its  maximum  value,  Gy  =  Ny  The  goal  is  to 
scan  the  superarray  focus  away  from  the  beacon  in  search  of  an  unknown 
signal  source  while  maintaining  an  acceptable  value  of  gain.  First 
each  subarray  scans  its  beam  to  the  location  S  as  shown  in  Fig.  6.1(c). 

To  focus  the  superarray  at  S,  the  phase  shift  determined  from  (5.6)  is 
applied  to  the  output  of  each  subarray  and  the  outputs  are  summed  as 
depicted  in  Fig.  6.1(d).  The  VLA  gain  at  S  is  determined  by  the  degrada¬ 
tion  of  coherence  due  to  the  random  fluctuations,  as  predicted  by  the 
MCF.  The  superarray  continues  to  scan  away  from  the  beacon  until  pair¬ 
wise  coherence  decreases  to  such  a  value  that  there  is  no  appreciable 
gain. 

Since  it  may  be  desirable  to  cover  a  larger  area,  it  is  necessary 
to  place  other  beacons  to  insure  continuous  coverage.  Each  beacon  has 
its  own  area  of  coverage,  and  the  beacon  locations  are  determined  by  the 
size  of  these  areas  so  that  coverage  js  continuous.  The  procedure  out¬ 
lined  above  is  then  repeated  for  each  beacon. 

It  is,  of  course,  necessary  that  the  required  density  of  beacons 
is  practical  for  the  given  system  specifications.  One  of  the  primary 
purposes  of  this  work  is  to  provide  a  procedure  for  determining  the 
feasibility  of  a  VLA  system  design  for  given  acoustic  parameters  and 
system  geometry,  within  the  limitations  of  the  oceanographic  fluctuations 
considered  here.  It  should  be  remembered  that  geographic  anomalies  have 
not  been  included  as  sources  of  fluctuations  and  will  be  a  source  of 
further  performance  degradation. 


AREA  OF  INCOHERENT  BEAM  INTERSECTION 


FULLY  COHERENT  FOCUS 


AREA  OF  INCOHERENT 
BEAM  INTERSECTION 


PARTIALLY  COHERENT  FOCUS 


Hg.  6.1  VLA  bcamfonn  and  scanning. 
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6.2  SYSTEM  DESIGN  PROCEDURE 


The  primary  considerations  in  the  design  of  a  VLA  of  subarrays 
are  the  performance  specifications  of  detection  ability  and  localiza¬ 
tion  ability.  Detection  ability  is  measured  by  the  system  gain  and 
is  determined  by  the  MCF.  Localization  ability  is  determined  primari¬ 
ly  by  the  system  configuration. 

The  primary  system  requirements  related  to  the  gain  are  the  num¬ 
ber  and  density  of  beacons  required  for  coverage  of  a  desired  area, 
given  the  system  configuration  and  the  acoustic  parameters  for  the 
ocean  area  of  interest.  Fig.  6.2  gives  the  value  of  the  MCF  required 
to  achieve  certain  values  of  VLA  gain,  G^,  as  a  function  of  the  num¬ 
ber  of  subarravs,  N„.  from  (2.63).  When  the  reauired  value  of  Y_  has 

V  *  £> 

been  determined,  the  area  of  coverage  with  one  beacon,  A^,  can  be 
found  from  the  contour  of  constant  coherence  using  the  results  of  Chap¬ 
ter  5.  It  was  shown  that  the  directions  of  extrema  of  scanning  ability 
are  approximately  parallel  and  perpendicular  to  the  VLA  baseline.  By 

computing  these  coherence  distances,  S  and  S  ,  respectively,  for  the 

x  y 

outermost  pair  of  receivers,  the  area  A  can  be  approximated  as  a 
rectangle. 


A 


B 


S  S 
x  y 


(6.1) 


Within  this  area  the  gain  will  exceed  the  minimum  required  value  since 
the  outermost  pair  of  receivers  has  the  lowest  coherence.  Assume  that 
each  beacon  in  the  area  of  interest  has  approximately  the  same  coher¬ 
ence  contour  with  the  same  area  A^.  Then  the  required  spacing  between 
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beacons  Is  S  In  the  direction  parallel  to  the  VLA  baseline,  and  S 
x  y 


In  the  perpendicular  direction.  For  a  desired  total  area  of  coverage, 
■A^,,  the  required  number  of  beacons  is 

(6.2) 


N 


B 


_T 

* 


Another  design  consideration  is  the  required  refocusing  time  for 
each  beacon.  In  Section  3.1  the  scanning  time  was  limited  to  T  «rpQ, 
where  Tq  is  a  characteristic  time  of  the  large  scale  environmental 
fluctuations,  which  is  on  the  order  of  days.  In  addition  there  was 
shown  to  be  no  dependence  on  scan  time  due  to  internal  waves  and 
tides.  For  internal  waves,  the  scan  time  determines  the  minimum  dis¬ 
tance  for  which  the  channels  are  independent;  thus,  for  scan  distances 
larger  than  this,  there  is  no  dependence  on  T.  Since  internal  tides 
were  shorn  to  have  a  negligible  effect  or  coherence,  their  dependence 
on  scan  time  can  be  ignored.  The  limiting  factor  on  scanning  time 


therefore  is  the  characteristic  time  T^. 


Assume  that  an  upper  limit,  Tg,  is  placed  on  scan  time  so  that 


Tg  «  rJIQ.  A  value  of  Tg  =  12  hr  may  be  reasonable,  but  due  to  the 


limited  knowledge  of  large  scale  fluctuations  it  should  be  determined 
by  experiment.  If  a  beacon  has  a  lifetime  T^,  and  if  >  Tg,  then 


each  subarray  must  refo  us  on  the  beacon  at  intervals  of  Tg.  However, 


if  Tn  <  T_  it  will  be  necessary  to  replace  each  beacon  at  intervals  of 
B  o 


T  .  This  is  an  important  consideration  for  system  design  implements- 
B 


tion  and  requires  further  study. 
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The  localization  ability  of  the  VLA  is  determined  primarily  by 
the  subarray  baamwidths  and  the  range  to  the  source.  If  each  subarray 
has  a, length  Lg, • then  the  beamwidth  is 


(6.3) 


When  the  separation  between  subarrays  is  large  the  area  of  intersection 
of  the  beams  at  a  range  R  is  then  approximately 


os  =  (KA8S)2 


(6.4) 


The  desired  resolution  determines  limits  on  the  relationships  between 
frequency,  range,  and  subarray  length.  The  number  of  resolution  cells 
per  beacon  is 


N 


R 


(6.5) 


A  requirement  for  feasibility  is  that  0g  <<  Ag  so  that  is  large. 

The  resolution  cell  of  the  VLA  focus  can  be  determined  from  (2.55) 


and  (2.57)  as 


ApBAS 


(6.6) 
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In  Fig.  6.3.  There  are  =  7  linear  subarrays  distributed  along  a 
baseline  of  Ly  =  150  km.  Each  subarray  has  Ng  ='40  sensors  spaced  one 
half  wavelength  apart  at  f  =  50  Hz  (X  =  30  m),  so  the  subarray  length 


is  Lg  =  585  ra.  If  the  noise  is  incoherent,  between  individual  sensors 
in  a  subarray,  then  the  subarroy  gain  is  0  =  16  dB  from  (2.62). 

It  is  desired  to  form  a  VLA  which  will  increase  the  system  gain 


by  a  minimum  of  =  6  dB  at  f  =  50  Hz.  The  desired  area  of  coverage 
2 

is  A^,  =  75000  km  centered  about  an  initial  beacon  range  of  Rq  =  250  km 
as  shorn  in  Fig.  6.4.  From  Fig.  6.2  the  required  value  of  the  KCF  is 
found  to  be  Yg  ~  0*5. 

Assume  that  the  multipath  parameters  are  K  =  16  and  0=2°.  To 

determine  the  scan  distances  S  and  S  ,  the  outer  pair  of  receivers 

x  y 

is  used  for  the  computation  since  they  vi]  ]  have  the  lowest  coherence. 

*  Based  on  calculations  using  random  array  theory. 


Hiis  Insures  that  the  gain  will  be  greater  than  6  dB  throughout  the  scan 

area,  A^.  At  and  near  the  beacon  the  gain  will  be  10log^Q7  =  8.5  dB. 

Thus  the  average  gain  within  the  y  »  0.5  contour  is  in  excess  of  7  dB 

and  the  maximum  gain  is  8.5  dB.  Using  tint  MCF  computation  procedure 

of  Section  5.3.2,  it  is  found  that  S  =  165  km  and  S  =  50  km,  giving 

x  y 

2 

a  total  coverage  area  of  An  =  8250  km  with  a  beacon  at  Rn.  Assuming 

B  u 

that  the  area  of  coverage  for  each  beacon  is  the  same,  the  total  number 
of  beacons  required  is  N  =  9,  from  (6.2).  The  beacon  configuration  and 
coverage  areas  are  illustrated  in  Fig.  6.5. 

The  subarray  beamwidth  is  found  to  be  A9g  =  .051  rad.  At 

2 

R  **  250  lan  the  resolution  cell  from  (6.4)  therefore  is  a  K  164  km  , 

U  D 

and  the  number  of  resolution  cells  per  beacon  is  N  =  50,  from  (6.5). 

K 

From  (6.6)  the  resolution  size  of  the  VLA  focus  is  found  to  be  Oy  - 
2 

4000  m  .  The  beacon  coverage  area  and  resolution  cell,  Og,  is  illus¬ 
trated  in  Fig.  6.6- 

In  summary,  this  VLA  will  increase  system  gain  by  more  than  6  dB, 

2 

covering  an  area  of  75000  km  with  9  beacons  spaced  by  50  km  in  the  per¬ 
pendicular  direction  and  165  km  in  the  parallel  direction.  The  size  of 

2 

the  resolution  cell  is  164  km  for  a  total  .of  50  resolution  cells  per 
beacon  and  450  resolution  cells  over  the  entire  coverage  area. 

6.3  CONSIDERATIONS  IN  SYSTEM  IMPLEMENTATION 

This  work  lias  been  primarily  concerned  with  the  most  important  VLA 
system  design  consideration,  that  of  signal  coherence  between  widely 
spaced  receivers.  The  derived  multipath  coherence  function  provides  the 


beacon  spacings  required  to  maintain  a  desired  coherence  and  VLA  gain. 
There  are  other  considerations  for  a  complete  system  design  implementa¬ 
tion  which  will  determine  system  feasibility.  A  detailed  discussion  of 
these  factors  is  oeyond  the  scope  of  this  work.  However,  this  section 
enumerates  the  most  important  of  them,  with  practical  suggestions  as  a 
basis  for  further  study. 

Beacon  placement 

For  a  practical  VLA  system  design  the  method  of  placement  of  bea¬ 
cons  is  an  important  consideration.  Permanent  beacon  installations 
would  be  expensive  with  a  lack  of  flexibility  in  location  and  a  high 
probability  of  discovery.  However,  temporary  beacons  with  a  limited 
lifetime  would  avoid  these  problems.  The  controlling  factors  in  deter¬ 
mining  the  feasibility  of  temporary  beacon  installations  would  be  the 
method  of  placement  and  the  beacon  lifetime,  T„.  Since  the  use  of 
temporary  beacons  implies  a  beacon  replacement  if  the  desired  scan  time 
about  a  beacon,  x„,  is  greater  than  its  lifetime,  T  ,  the  method  of 

b  IS 

beacon  placement  should  be  expedient.  One  method  that  warrants  consider¬ 
ation  is  dropping  beacons  from  an  aircraft.  The  technology  in  this 
area  is  well  developed  and  the  method  offers  the  obvious  advantages  of 
flexibility  in  beacon  location,  acciu..rv  of  location  by  navigational 
methods,  and  ease  of  immediate  beacon  replacement. 

Another  possibility  is  the  use  of  beacons  of  opportunity  such  as 
surface  shipping  [1].  The  advantages  arc  availability  at  no  expense, 
and  concealment.  There  will  be  difficulties  in  phase  measurement  due  to 
spatial  variations,  and  a  lack  of  reliability  and  flexibility  in  location. 
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However,  due  to  the  abundance  of  shipping  traffic  in  areas  where  hoar- 
con  placement  would  be  difficult  or  impossible,  it  should  be  explored 
as  a  possibility. 

Beacon  waveforms 

The  most  important  requirement  for  a  beacon  waveform  is  that  it 
enables  accurate  measurement  of  the  channel  impulse  response  in  the 
presence  of  noise.  The  waveforms  of  all  beacons  must  be  known  at  all 
receiver  sites  and  they  must  be  distinguishable.  In  addition,  the  bea¬ 
con  signals  should  be  undetectable  to  all  others. 

Measurements  of  the  impulse  response  of  underwater  channels  and 
acoustic  phase  detection  have  been  investigated  theoretically  [2], 

[3],  [4]  and  experimentally  [51  -  [10].  It  appears  that  pseudo¬ 
noise  (PN)  sequences  [11]  satisfy  the  requirements  stated  above  and 
should  be  considered  in  a  superarray  system  design. 

Some  of  the  characteristics  which  the  PN  code  should  possess  will 
be  dictated  by  the  channel  characteristics.  The  time  length  of  the 
code  must  be  greater  than  the  time  spread  of  the  channel  to  insure 
unambiguous  measurement  of  the  multipath  arrivals.  The  sequences  must 
be  distinguishable  between  beacons,  thus  a  different  code  should  be 
used  for  each  beacon.  Each  pair  of  sequences  should  have  good  cross¬ 
correlation  properties  so  that,  only  the  desired  beacon  waveform  is  de¬ 
tected.  Another  consideration  is  the  time  required  for  each  receiver 
to  synchronize  with  the  beacon  PN  sequence. 


Source  localization 
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The  VLA  localization  accuracy  is  primarily  determined  by  the 
accuracy  of  location  of  the  beacons  and  subarrays,  and  by  the  number 
of  ambiguous  VLA  focal  areas  within  the  area  of  intersection  of  the 
subarray  beams.  If  there  are  VLA  ambiguities  within  cf  ,  then  Cc  is 
the  minimum  resolution  cell.  The  linear  dimension  of  ag  is  typically 
on  the  order  of  tens  of  kilometers.  The  beacon  locations  and  the  loca¬ 
tions  and  orientations  of  the  subarrays  will  be  known  to  at  least 
navigational  accuracy,  whose  error  is  much  less  than  this  value. 
Therefore,  it  can  be  assumed  that  there  is  little  effect  on  the  local¬ 
ization  accuracy  determined  by  the  resolution  cell  a g. 

It  was  shown  in  the  previous  section  that  if  there  are  no  VLA 

ambiguities  within  crg,  then  the  size  of  the  resolution  cell  will  be 

2 

decreased  to  which  is  on  the  order  of  A  .  However,  the  location 
of  is  highly  sensitive  to  the  location  accuracy  of  the  system  com¬ 
ponents.  Even  if  the  subarray  and  beacon  locations  could  be  known 
within  fractions  of  a  wavelength,  the  location  of  would  still  be  in. 
error  due  to  the  randomness  of  the  medium.  However,  due  to  the  poten¬ 
tial  increase  by  orders  of  magnitude  in  localization  ability,  this  sub¬ 
ject  should  receive  further  study. 

Source  motion 

An  application  of  a  VLA  system  to  the  detection  of  moving  sources 
presents  additional  complications.  In  Section  ?. .4  it  was  shown  that 
travel  time  fluctuations  must  vary  slowly  compared  to  signal  duration 
time,  so  that  the  channel  transfer  functions  will  be  time  invariant. 
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However,  this  may  not  be  valid  for  a  moving  source,  because  the  spatial 
fluctuations  due  to  multipath  interference  will  vary  with  time.  The 
seriousness  of  this  effect  will  depend  upon  integration  time  and  source 
velocity.  A  related  problem  is  the  ability  to  track  the  source  by  main¬ 
taining  the  VLA  focus  on  its  changing  location.  The  source  motion  also 
causes  a  complicated  ^oppler  effect  due  to  a  different  frequency  shift 
in  each  ray.  However,  this  effect  can  be  minimized  by  properly  shifting 
the  center  frequency  of  the  receiver  filters. 

Post-detection  focusing  and  tracking 

After  initial  detection  of  a  signal  source  with  tie  VLA,  it  may  be 
possible  to  further  increase  the  signal  to  noise  ratio  by  enhancing  e 
partially  coherent  VLA  focus.  Each  subarray  would  measure  the  relc._j.ve 
signal  phase  or  coherence  of  the  signal  waveforms.  Using  this  informa¬ 
tion  a  refocused,  high  resolution  spot  is  placed  on  the  source  by  self- 
cohering  or  adaptive  bearaforming  techniques.  The  focus  is  then  scanned 
in  the  vicinity  of  the  source  for  the  purpose  of  tracking. 

Geographic  fluctuations 

In  Chapter  3  geographic  fluctuations  such  as  currents  and  eddies 
were  discussed.  The  theoretical  developments  in  this  work  were  limited 
to  those  which  are  not  geographic  in  nature.  However,  due  to  the 
prevalence  of  geographic  anomalies,  they  must  be  considered  in  a  VLA 
system  design. 

Due  to  the  variability  and  unpredictability  of  some  of  these  fluc¬ 
tuations  it  is  difficult  to  determine  their  effect  on  coherence.  The 
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best  approach  would  be  to  evaluate  the  effect  of  geographic  fluctua¬ 
tions  by  experiment  in  the  ocean  area  of  interest  for  a  VLA  system. 

The  geographic  areas  of  prevalence  of  some  of  these  fluctuations  might 
be  considered  in  the  geographic  location  of  a  VLA. 

Determination  of  actual  phase  for  VLA  scanning 

The  development  of  the  MCF  in  Chapter  4  theoretically  predicted 
the  average  phase  shift  required  for  each  subarray  in  order  to  scan 
the  VLA  with  partial  coherence,  however,  this  result  depends  on  an 
accurate  knowledge  of  the  multipath  structure  which  may  not  be  available. 
It  is  important  to  know  the  correct  average  phase  for  each  individual 
situation.  If  the  average  phase  shift  is  inaccurate  then  another  ran¬ 
dom  variable  i?  Introduced  which  will  further  degrade  coherence. 

This  suggests  the  desirability  of  experimentally  measuring  phase 
as  a  function  of  scan  distance.  This  measurement  will  show  a  phase 
trend  [12]  with  fluctuations  about  the  trend  due  to  the  variations  con¬ 
sidered  in  this  work.  This  procedure  of  surveying  the  scan  area  is 
performed  only  once,  and  the  phase  trend  measured  is  then  used  as  the 
average  phase  shift  for  future  VLA  scanning.  The  true  phase  will  vary 
causing  a  degradation  of  coherence,  but  the  trend  should  remain  constant. 

Coherent  noise  sources 

In  Section  2.3  noise  was  limited  to  random  broadband  ambient  noise 
which  is  incoherent  between  VLA  subarrays.  However  there  is  a  possi¬ 
bility  of  discrete  shipping  interference  which  may  be  coherent  between 
subarrays.  Experimental  [13]  and  theoretical  [14]  results  can  be  used 


to  predict  this  shipping  density  for  the  North  Atlantic.  A  method  of 
near  field  adaptive  nulling  of  coherent  noise  with  a  VLA  of  subarrays 
was  developed  in  [14]  based  upon  the  concept  of  null  steering  [15]. 

Practical  implementation  of  this  technique  would  involve  an  initial 
localization  of  interfering  shipping  by  airborne  radar  detection  or 
other  means,  and  a  null  tracking  system  in  each  subarray  so  that  individ¬ 
ual  nulls  in  the  subarray  patterns  can  follow  the  shipping  traffic. 

The  near  field  pattern  of  the  VLA  can  then  be  visualized  as  having 
"holes”  which  follow  the  ships  as  they  move  throughout  the  area. 

Subarray  location 

Some  additional  system  flexibility  can  be  acquired  if  the  subarrays 
can  be  placed  in  arbitrary  locations.  A  possible  VLA  system  might  con¬ 
sist  of  several  floating  random  arrays  [16]  which  could  be  deployed  by 

* 

an  aircraft  in  any  desired  locations  .  Combined  with  the  use  of  air¬ 
dropped  beacons,  the  VLA  system  would  then  have  the  advantage  of  com¬ 
plete  mobility.  The  disadvantage  would  be  a  further  degradation  of  sys¬ 
tem  gain  due  to  the  larger  spacings  between  sensors  in  a  random  float¬ 
ing  array. 

6.4  SUMMARY 

This  chapter  has  presented  an  application  of  the  MCF  to  the  design 
of  a  VLA  of  widely  separated  subarrays.  The  general  system  design 

*This  idea  was  suggested  by  Professor  F.  Labor,  Moore  School  of  Electrical 
Engineering,  University  of  Pennsylvania. 


approach  was  outlined.  A  procedure  was  then  developed  for  determining 
beacon  spacings  required  for  a  given  VIA  configuration  to  maintain  a 
specified  gain  over  a  desired  coverage  area.  VLA  refocusing  times  were 
shown  to  be  dependent  on  the  large  scale  oceanographic  fluctuations. 
Localization  ability  was  discussed  in  terms  of  the  subarray  beamwidths 
and  the  size  of  the  focal  area  of  the  VLA. 

A  design  example  was  presented  for  some  realistic  system  para¬ 
meters.  This  example  showed  that  coherent  combination  of  7  subarrays 

could  increase  system  gain  by  an  average  of  more  than  7  dB  over  an  area 
2 

of  75000  km  with  the  use  of  only  9  beacons. 

Finally,  some  important  considerations  in  system  implementation 
were  discussed,  and  proposals  were  made  for  practical  solutions. 
Specifically  mentioned  were  the  possibilities  of  beacon  placement  by 
aircraft,  and  PN  sequences  for  beacon  waveforms.  The  idea  of  floating 
subarrays  also  deployed  from  an  aircraft  was  discussed  as  a  method  of 
making  an  entire  VLA  mobile. 
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CHAPTER  7 


SUMMARY  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

7.1  SUMMARY  OF  RESULTS 

In  order  to  view  the  results  of  this  work  in  the  proper  per¬ 
spective,  it  is  helpful  to  review  the  line  of  reasoning  that  led  to 
their  development.  The  motivation  for  this  work  was  the  idea  of 
coherently  combining  widely  spaced  subarrays  in  a  random  multipath 
underwater  medium.  The  purpose  of  forming  this  very  large  array  is 
to  increase  the  potential  signal  to  noise  ratio  and  the  localization 
ability.  The  enhancement  of  detection  ability  is  measured  by  the 
array  gain,  defined  in  terms  of  the  signal  coherences  between  all 
pairs  of  subarrays.  The  foremost  problem,  then,  was  to  develop  a 
solution  for  this  coherence  in  terms  of  the  environmental  and  acoustic 
characteristics  of  the  ocean. 

This  led  to  the  development  of  a  new  measure  of  coherence,  called 
the  multipath  coherence  function,  defined  in  terms  of  ensemble  averages 
of  the  random  transfer  functions  of  the  multipath  channels.  Since  the 
receivers  are  widely  spaced  the  channels  are  stochastically  indepen¬ 
dent.  This  important  simplification  demonstrated  the  existence  of 
coherence  without  correlation  between  channels;  it  also  enabled  the 
MCF  to  be  factored  into  separate  auto-coherences,  making  the  final  solu 
tion  mathematically  feasible.  Another  important  simplification  was  the 


envelope  approximation  for  the  auto-coherence,  which  factored  the 
effects  of  random  fluctuations  from  those  of  frequency  selective  multi- 
path  interference.  The  MCF  was  then  formulated  as  a  function  of  source 
range  and  scanning  distance,  for  general  oceanographic  fluctuations. 

It  then  remained  to  specify  the  stochastic  parameters  of  the  MCF 
for  real  oceanographic  fluctuations.  This  required  original  analyses 
of  the  effects  of  spatial  variations  due  to  multipath  interference,  and 
of  internal  tides,  on  coherence.  The  stochastic  parameters  of  internal 
wave  fluctuations  were  obtained  from  the  literature.  A  comparison  of 
these  effects  then  showed  that  spatial  variations  were  predominant  in 
scanning,  while  internal  tides  were  a  negligible  influence. 

The  remaining  step  was  to  apply  these  results  to  the  initial  objec¬ 
tive  of  predicting  VLA  performance  in  terms  of  signal  to  noise  gain  for 
given  system  configurations.  The  system  design  approach  was  to  use 
self-cohering  techniques  whereby  the  VLA  initially  focuses  on  a  known 
beacon  source  in  the  near  field,  and  then  scans  in  the  vicinity  of  the 
beacon  in  search  of  an  unknown  signal.  Thus  the  quantities  of  interest 
were  the  number  and  spacing  of  beacons  required  to  maintain  a  specified 
gain  while  scanning  the  VLA  between  beacons,  A  design  example  for  some 
realistic  parameters  then  showed  the  existence  of  significant  coherence 
over  large  ocean  areas.  The  conclusion  is  that  a  VLA  design  might  be 
possible  and  practical. 

In  summary,  there  are  three  primary  results  from  this  research. 

The  first  is  a  general  solution  for  signal  coherence  in  uncorrclated 
multipath  channels.  The  second  is  a  specific  application  of  this 
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multipath  coherence  function  to  the  design  of  a  VLA  composed  of  wide¬ 
ly  spaced  subarrays.  Finally ,  numerical  .results  showed  that  such  a 
VLA  design  is  feasible  for  certain  system  configurations  and  multipath 
characteristics. 


7.2  ANALYSIS  OF  RESULTS 

The  following  sections  outline  the  important  conclusions  to  be  made 
from  the  results  of  this  work,  and  an  explanation  of  its  limitations. 

The  points  considered  are  limited  to  the  areas  of  the  three  primary 
results  of  this  work  stated  in  Section  7.1.  Further  information  and  de¬ 
tail  can  be  obtained  from  the  summaries  at  the  end  of  relevant  chapters. 


7.2.1  CONCLUSIONS 

The  most  important  conclusion  to  be  made  from  the  development  of  the 
Col^A-e^C-C  FwcSYo* 

_ _ — £MC^)  is  that  it  demonstrates  the  existence  of  coherence  without  correla¬ 
tion  between  random  channels.  The  MCF  demonstrates  the  importance  of 
the  size  of  fluctuations  compared  to  their  correlation.  The  existence 
of  partial  coherence  implies  a  non-zero  mean  signal  field  for  fluctua¬ 
tions  which  are  small  enough.  It  was  also  shown  that  the  MCF  is  indepen¬ 
dent  of  the  signal  source  and  depends  only  on  the  properties  of  the 

/ 

medium. 

The  importance  of  frequency  dotaain  processing  is  readily  observed 
by  comparing  the  coherence  function  with  the  normalized  cross-correlation 
function.  A  broad  band  signal  waveform  in  a  multipath  medium  may  have 


only  one  or  two  discrete  frequencies  at  which  coherence  is  high.  The 
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cross-correlation  function,  however,  considers  the  entire  signal  wave¬ 
form  and  will  have  a  much  smaller  value  than  the  maximum  value  of  the 
coherence  function. 

-  The  relevance  of  the  MCF  is  in  its  relation  to  array  gain.  The 
,  magnitude  of  the  MCF  gives  the  signal  power  gain  achieved  by  combining 

I  the  outputs  of  a  pair  of  receivers.  Its  phase  is  the  average  phase 

i 

I  difference  between  the  received  signals  which  is  required  to  combine 

I  -  ' 

them  with  partial  coherence. 

The  mathematical  solution  for  coherence  led  to  a  convenient 
factorization  into  nine  auto-coherence  terms.  The  stochastic  indepen¬ 
dence  of  the  receiver  channels  permits  an  auto-coherence  to  be  com¬ 
puted  for  each  channel  independent  of  the  others  resulting  in  a  large 
computational  savings.  The  envelope  approximation  further  factored 
the  effects  of  random  fluctuations  from  those  of  frequency  selective 
multipath  interference.  This  allows  the  prediction  of  maximum  coherent 
frequencies  independent  of  the  actual  multipath  ray  configuration.  The 
extension  of  the  MCF  to  VLA  scanning  led  to  a  further  factorization  of 
the  MCF  into  an  auto-coherence  due  to  the  effects  of  spatial  variations 
in  the  multipath  interference.  The  advantages  of  these  factorizations 
are  computational  simplicity  and  the  ability  to  compare  various  effects 
on  coherence  readily. 


The  MCF  formulation  is  a  simple,  concise  mathematical  expression, 
and  does  not  depend  on  present  knowledge  of  oceanographic  fluctuations. 
The  solution  is  adaptable  to  future  developments  in  the  causes  of  these 
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fluctuations. 
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The  complete  value  of  the  MCF  is  determined  by  multiplying  the 
eight  auto-coherence  factors  for  a  given  channel,  system  geometry,  and 
acoustic  frequency.  It  has  been  shown  that,  for  a  constant  value  of 
O,  and  for  a  coherent  frequency  (y  =  1),  the  MCF  increases  as  the  num¬ 
ber  of  rays,  K,  increases.  However,  if  an  increase  in  the  number  of  ray 
arrivals  causes  a  corresponding  increase  in  the  ray  spread,  a,  the  coher¬ 
ence  may  decrease.  This  fact  is  an  important  consideration  in  the 
choice  of  receiver  location.  For  example,  it  may  be  wise  to  place  the 
receivers  at  depths  where  there  is  a  large  number  of  ray  arrivals  with¬ 
in  a  small  angular  spread,  rather  than  to  choose  a  location  having  only 
one  ray  with  the  hope  of  avoiding  spatial  variations  due  to  multipath 
interference  entirely. 

The  identification  of  the  auto-coherence  with  specific  oceanographic 

fluctuations  allowed  the  relative  effect  of  each  type  of  fluctuation  tc 

be  determined.  For  each  receiver  there  is  an  auto-coherence  due  to 

internal  waves  in  the  beacon  channel,  ytI;  an  auto-coherence  due  to 

w 

internal  waves  and  spatial  multipath  interference  in  the  scan  channel, 

Yy  g5  and  an  auto-coherence  due  to  frequency  selective  multipath  inter¬ 
ference,  Yjj.  There  is  also  the  effect  of  internal  tides  on  the  coherence 
between  receivers,  y^. 

The  factor  y^  depends  upon  the  typical  angle  at  which  rays  cross 
the  sound  channel  axis  and  has  a  higher  value  for  steeper  rays.  It 
decreases  with  the  source  range,  R,  and  with  the  acoustic  frequency,  f. 

Y'j  g  contains  the  effect  of  internal  waves  in  the  scan  channel  plus 
the  effect  of  spatial  variations  due  to  multipath  interference  in 
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scanning.  The  effect  of  spatial  variations  depends  upon  the  change 
in  source  range,  x,  in  scanning  away  from  the  beacon.  Thus,  the  effect 
is  most  severe  when  scanning  in  a  direction  perpendicular  to  the  base¬ 
line  of  a  VLA.  The  parameter  which  affects  the  spatial  variations  is 
the  angular  ray  spread  at  the  receiver,  a,  and  it  was  found  that  coher¬ 
ence  decreases  as  <7  increases.  For  a  given  value  of  O  the  coherence 
decreases  in  scanning  as  a  function  of  |x|/a. 

The  effect. on  coherence  due  to  frequency  selective  multipath 
interference,  y.  ,  was  found  to  depend  on  the  nominal  time  configura¬ 
tions  of  the  multipath  arrivals  for  each  receiver  channel.  The  coher¬ 
ence  Yjj  is  specified  primarily  in  terms  of  coherent  frequencies,  i.e., 
frequencies  at  which  the  rays  interfere  constructively,  and  the  coherent 
bandwiciths  cantered  on  these  frequencies.  For  a  simplified  equal  time 
spacing  formulation  of  the  ray  arrivals,  it  was  found  that  coherent 
frequencies  occur  at  harmonics.  The  fundamental  coherent  frequency  in¬ 
creases  as  the  number  of  rays,  K,  increases,  and  as  the  time  spread, 

T_,  decreases.  The  coherence  bandwidths  are  inversely  proportional  to 

u 


It  was  found  that  internal  tides  had  a  negligible  effect  on  coher¬ 
ence  in  scanning  compared  to  the  effect  of  spatial  multipath  interfer¬ 
ence.  Consequently  y^,  has  a  value  of  unity  for  all  scan  distances  of 
practical  interest.  However,  the  coherence  due  to  frequency  selective 
multipath  interference  is  a] so  a  consideration,  since  the  spacing  between 
coherent  frequencies  increases  as  K  increases.  All  of  these  factors 
should  be  taken  into  account  for  an  optimum  system  design. 
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The  argument  of  the  MCF  is  the  average  phase  difference  between 
signals  necessary  to  combine  them  with  partial  coherence.  The  only 
.contribution  to  this  phase  is  due  to  the  spatial  multipath  interference. 
It  was  found  that  the  phase  has  an  expected  linear  variation,  k^x,  and 
a  contribution  due  to  the  ray  spread,  a.  It  was  found  that  the  devia¬ 
tion,  of  the  average  phase  from  the  linear  trend  is  less  than  n/4  radians. 

The  actual  design  of  a  superarray  system  requires  the  use  of  bea¬ 
cons  with  known  locations  and  known  waveforms  for  the  purpose  of 
initially  focusing  the  VLA  due  to  the  unknown  state  of  the  medium. 
Scanning  th,'  VLA  is  performed  by  first  scanning  the  subarray  beams  to 
the  desired  location  and  then  applying  the  required  phase  shifts  to  the 
subarray  outputs.  These  phase  shifts  are  nominal  or  average  values 
which  are  either  predicted  from  a  ray  tracing  program  or  experimentally 
measured.  The  MCF  then  predicts  the  defoc using  of  the  VLA  due  to  the 
fluctuations  about  the  average  phase.  The  localization  accuracy  of  the 
VLA  is  determined  primarily  by  the  area  of  intersection  of  the  subarray 
beams  because  of  ambiguities  in  the  VLA  pattern. 

The  MCF  predicts  the  contours  of  constant  coherence  giving  the 
area  of  coverage  with  one  beacon  for  a  desired  array  gain.  This  deter¬ 
mines  the  number  of  beacons  and  their  geometric  configuration  for  a  re¬ 
quired  total  area  of  coverage. 

Numerical  results  for  a  specific  multipath  configuration  and  a  VLA 

of  7  subarrays  predicted  an  average  gain  in  excess  of  6  dB  over  an  area 
2 

of  75000  W  with  the  use  of  only  9  beacons. 


7.2.2  LIMITATIONS 


The  restrictions  of  the  foregoing  theory  are  few*  and  are  simply 
stated.  There  are  some  limitations,  however,  which  may  be  important 
considerations  in  a  VLA  system  design,  and  are  itemized  below. 

•  Geometrical  optics  -  the  theory  of  fluctuations  has  been  limited 
to  the  geometrical  optics  regime  (small  fluctuations),  with 
associated  limits  on  frequency  and  range.  However,  since  the 
fluctuations  are  uncorrelated,  the  MCI'  depends  only  on  their  size. 
Consequently,  coherence  would  be  very  small  for  fluctuations  larger 
than  those  of  the  geometrical  optics  region,  so  it  is  unnecessary 
to  consider  the  other  regions  of  fluctuations. 

•  Small  scale  size  of  fluctuations  -  the  limitation  to  fluctuations 
of  small  correlation  distance  and  time  (internal  waves)  places  re¬ 
strictions  on  scan  distance  and  scan  time.  But  since  coherence 

is  low  for  larger  scan  distances  due  to  the  fluctuations  considered 
here,  there  is  no  need  to  consider  larger  scale  size  fluctuations. 

•  Non-geo graphic  fluctuations  -  the  theory  here  does  not  consider 
geographic  anomalies  such  as  currents  and  eddies.  These  fluctua¬ 
tions  should  be  seriously  considered  in  a  VLA  system  design,  either 
by  measuring  their  effect  for  the  area  of  interest,  or  by  avoiding 
them  entirely  in  locating  the  VLA. 

•  Source  motion  -  a  restriction  imposed  early  in  this  work  was  that 
signal  duration  tine  be  much  less  than  the  characteristic  time  of 
all  fluctuations.  Source  motion  implies  the  existence  of  spatial 
v..  iatious  of  multipath  interference  in  the  signal  duration  time. 
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There  is  therefore  a  requirement  that  the  source  motion  is  minimal- 
for  the  time  that  it  takes  the  signal  processor  to  compute  the 
power  spectral  densities.  The  effect  of  Doppler  can  be  accounted 
for  by  simply  shifting  the  center  frequency  of  the  signal  processor 
filter.  This  subject  requires  further  study. 

Horizontal  scanning  -  although  no  restriction  was  placed  on  receiver 
location  due  to  the  use  of  a  beacon  for  initial  cophasing,  scanning 
was  limited  to  a  horizontal  plane  at  the  beacon  depth.  It  is 
assumed  that  there  are  no  great  variations  in  depth  for  typical 
signal  sources.  The  vertical  coherence  distance  therefore  can  be 
presumed  to  be  large  enough  to  detect  all  sources  of  interest  with 
beacons  at  only  one  or  two  depths. 

Accuracy  of  beacon  and  receiver  locations  -  the  VI, A  localization 
ability  is  determined  by  the  area  of  intersection  of  the  subarray 
beams,  which  may  be  on  the  order  of  tens  of  square  kilometers. 

Since  receiver  and  beacon  locations  will  be  known  within  areas  much 
less  than  this,  there  will  be  essentially  no  effect  on  localization 
ability  of  the  VLA. 

Incoherent  noise  -  if  the  noise  is  incoherent  between  receivers 
there  will  be  no  effect  on  VLA  gain  due  to  ambiguities  in  the  re¬ 
ceiving  pattern.  Coherent  noise  sources,  such  as  shipping  traffic, 
will  cause  a  decrease  in  gain  if  located  at  one  of  the  ambiguous 
focal  points  of  the  VLA.  The  procedure  of  near  field  adaptive 
nulling  discussed  in  Chapter  6  may  avoid  this  problem. 


•  Isotropic  multipath  field  -  it  was  assumed  that  the  number  of 
rays,  the  nominal  travel  times,  and  the  angles  of  arrival,  are 
invariant  throughout  the  scan  area  for  each  beacon.  The  source 
ranges  and  scan  areas  for  which  this  condition  is  fulfilled  must 
be  predicted  from  a  ray  tracing  program.  This  is  a  further  consid¬ 
eration  in  beacon  and  receiver  location. 

7.3  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Section  6.3  has  discussed  other  considerations  in  a  VLA  system  design 
implementation,  and  Section  7.2.2  has  stated  limitations  of  the  theory 
and  its  applications.  The  topics  mentioned  warrant  further  study  and  are 
summarized  here  vith  further  recommendations  concerning  experimental 
verification  of  results. 

1  -  An  experimental  test  of  the  theory  of  the  MCF  using  beacons. 

2  -  An  experimental  test  of  application  to  a  VLA  system  design. 

3  -  A  study  of  beacon  waveforms  such  as  PN  sequences. 

4  -  Evaluation  of  methods  of  beacon  placement  such  as  beacons  of 

opportunity  and  air-dropped  beacons. 

5  -  An  experimental  test  of  source  localization  ability. 

6  -  A  study  of  the  effects  of  source  motion. 

7  -  A  study  of  the  effects  of  geographic  fluctuations  such  as 

currents  and  eddies. 

8  -  Evaluate  the  effects  of  coherent  noise  sources. 

9  -  A  study  of  post-detection  focusing  and  tracking. 


Consideration  of  methods  of  subarray  location  such  as  the 

i 

use  of  floating  random  arrays. 

Study  of  the  use  of  sources  of  opportunity  instead  of 
beacons. 


APPENDIX 


THE  ENVELOPE  APPROXIMATION 

It  Is  desired  to  approximate  the  envelope  of  the  function 


y2(w)  = 


Kc  (u)  |  IIq  (w)  | 


l+[KjlIn((j)j  -l]c  (o) 


(A.l) 
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In  the  above  equation  ]Hq(o))|  is  the  square  of  the  normalized  multi- 
path  transfer  function  exhibiting  frequency  selective  fading.  It  is 
shown  in  Section  5.2.4  that  coherent  frequencies  located  at  the  primary 
maxima  of  the  pattern  have  typical  spacings  on  the  order  of  1  Hz,  and 


that  the  nulls  adjacent  to  each  coherent  frequency  are  separated  by 


fractions  of  1  Hz.  Between  coherent  frequencies  there  are  a  number  of 


secondary  maxima  which  have  generally  the  same  spacings  between  their 

9 


adjacent  nulls.  The  function  c“(w)  is  the  squared  characteristic  func¬ 
tion  of  random  fluctuations  and  is  a  monotcnically  decreasing  function 
of  frequency.  The  typical  characteristic  functions  considered,  that  of 


internal  vaves  in  (5.3)  and  of  spatial  variations  in  (5.6),  vary  slow- 

1 2 


ly  compared  to  |Hq(w)|  ,  and  can  be  considered  constant  between  any  two 
nulls. 

2  2  2 


At  the  nulls  of  ]Hq(w)J  ',  y  (w)  =  0  when  c  (w)  <  1,  so  the  local 


maxima  of  y  (w)  also  occur  between  the  nulls.  Consider  a  primary  mnxi- 

1 2  ^  ..  __  ..  _  t ..  /  .  v  i  2 


mum  of  |Hq(u>)|  located  at  w  ~  Wq  so  |Hq(u>q)|“  =  1.  Then  in  the  region 
between  the  adjacent  nulls 


tr  , 


158 


i 


159 


Y2(w)  = 


Kc2(a>0)lHQ(uj)l2 

l+lKjH0(w)|2-l]c2(w0) 


(A.  2) 


2  2 

Since  y2(w)  increases  monotonically  with  |Hq(u)]  when  c  (Wq)  <  1, 
its  maximum  value  in  the  region  between  the  nulls  about  wQ  is 


max 


V2(w) 


Kc2 (w0) 
l+(K-l)c2  ((jJq) 


(A. 3) 


2 

This  is  true  for  each  primary  maximum,  so  the  envelope  of  Y  (u)  f°r  a-^ 
u)  is 


onv  Y2  =  '4  " 


Kc2  (w) 


l+(K-l)c  (ca) 


(A.  A) 


A 

Since  the  primary  maxima  and  zeroes  of  Y  (w)  coincide  with  those  of 
o  2 

|Hq(w)|  the  approximation  to  Y  (w)  Is  now  written  as 

Y2(w)  =  y|(w)Y2j(w)  (A*5 


where 


y2j(w)  "  ly«>l2  • 


(A.  6) 


The  fractional  error  in  this  approximation  is 


e  " 


(A. 7) 


Equation  (A.l)  can  be  rewritten  as 


Y 


2 


1 

l-U-T&Y 


(A. 8) 


Then 


e 


(I-Y^Ye 


(A.  9) 


2 

At  the  primary  maxima,  y^  =  1,  so  there  is  no  error.  These  are  the 

regions  of  main  interest  since  they  are  the  locations  of  coherent 

frequencies.  The  error  in  the  approximation  is  greater  at  secondary 

2  1 

maxima.  For  example,  consider  a  secondary  maximum  where  y  =  y  . 

Assume  that  y2  =  -r  .  Then  y2  =  y,  but  y2  =  y,  so  £  =  25%.  The 

2 

fractional  error  continues  to  increase  as  y^  decreases. 

It  should  be  noted  that  the  envelope  approximation  is  not  valid 

2 

when  there  are  no  random  fluctuations,  i.e.,  when  c  (w)  ■=  1  for  all  w 

2 

since,  from  (A.l),  y'(w)  =  1  for  all  w. 
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